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PREFACE 
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SUMMARY 


The  primary  purpose  of  this  research  was  to  measure  experimentally  the 
initial  amplitude  and  growth  of  instability  waves  In  a  laminar  water  boundary 
layer  under  the  influence  of  an  active  wall.  The  research  was  conducted  in 
support  of  the  Compliant  Coating  Drag  Reduction  (CCDR)  Program  for  the  Office 
of  Naval  Research  (ONR).  Although  the  final  objective  was  not  met  as  a  conse¬ 
quence  of  limitation  in  funds,  several  goals  were  achieved  in  its  support.  The 
publications,  reports,  and  presentations  for  this  program  are  listed  in  the 
following  sections,  and  the  technical  accomplishments  are  described  in  detail 
in  the  appendices  which  include  all  publications  and  reports. 

Appendix  D  provides  the  theoretical  basis  for  the  feasibility  of  an 
active-wall  device.  The  analysis  involves  a  Kelvin-Helmholtz  solution  for  flow 
over  a  semi-infinite  surface  with  traveling  surface  waves.  Although  the  leading 
edge  of  the  surface  generates  a  secondary  flow,  wall  motion  can  be  adjusted  to 
modify  or  eliminate  freestream  disturbances  and  the  secondary  flow. 

Appendix  A  describes  a  noncontacting  electrooptic  displacement  sensor  for 
the  measurement  of  active-wall  motion.  The  sensor  monitored  surface  displace¬ 
ments  at  a  stand-off  distance  of  23  cm  with  a  two-element  position  sensitive 
detector  and  a  two-arm  optical  tri angulation  method.  The  displacement  sensor, 
which  was  designed  for  operation  in  either  air  or  water,  demonstrated  a  resolution 
of  better  than  1.3  pm  over  a  range  of  1.3  mm.  The  system  was  tested  with  two  light 
sources,  an  infrared  LED  (light  emitting  diode)  and  a  helium-neon  laser. 

Appendix  B  presents  the  results  on  the  development  of  an  active  wall 
which  was  an  all  solid-state  electromechanical  device  for  the  generation  of  small 
traveling  surface  waves.  The  active  wall  was  nominally  designed  for  the  cancel¬ 
lation  of  Tollmien-Schlichting  waves  in  a  laminar  boundary  layer  in  water.  The 
active  surface  area  of  the  prototype  device  was  51  mm  wide  by  40  mm  long  in  the 
flow  direction.  Wave  form,  wave  speed,  frequency,  and  amplitude  of  the  mechanism 
were  independently  selectable.  The  electronic  subsystem  was  a  hybrid  analog/digital 
device,  and  the  wall  motion  was  produced  by  piezoelectric  ceramic  elements.  Sinu¬ 
soidal  wall  motion  was  measured  in  air  over  an  amplitude  range  of  0.6  -  13  pm  and 


frequency  span  of  10  -  150  Hz.  The  maximum  wave  speed  was  166  cm/s.  The  paper 
also  includes  a  review  of  previous  active-wall  devices. 

Finally,  Appendix  C  presents  the  details  on  the  design  of  a  low- turbulence 
water  tunnel  for  an  active  wall  experiment.  Although  the  test  section  was  de¬ 
signed  for  an  active-wall  experiment,  other  laminar  and  turbulent  boundary  layer 
experiments  can  be  performed  in  the  facility.  The  water  tunnel  was  designed  for 
a  velocity  of  10  m/s  (32.8  ft/s)  in  the  test  section  whose  height,  width,  and 
length  are,  respectively,  15.2,  22.9,  and  76.2  cm  (  6  x  9  x  30  in).  The  top  and 
two  sides  of  the  test  section  were  constructed  from  Plexiglas  for  flow  visualiza¬ 
tion  and  laser  Doppler  (LDA)  measurements  while  the  bottom  surface  is  a  smooth 
flat  aluminum  plate  for  boundary  layer  measurements.  Estimates  indicate  that 
the  longitudinal  relative  turbulence  intensity  in  the  test  section  should  be  less 
than  0.2%.  Although  the  tunnel  has  been  assembled,  the  flow  characteristics  of 
the  test  section  have  not  been  measured. 

The  water  tunnel  was  fabricated  and  assembled  with  funds  from  Southwest 
Research  Institute  (SwRI).  Most  of  the  water  tunnel  was  constructed  from  stain¬ 
less  steel.  The  contraction  section  is  a  matched-cubic  curve  with  a  polished 
stainless  steel  interior,  and  its  contraction  ratio  is  16.7.  The  settling  chamber 
contains  six  fine-mesh  screens  and  stainless  steel  honeycomb.  The  facility  also 
includes  an  in-line  aluminum  heat  exchanger  which  also  functions  as  a  flow 
straightner. 
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A  NONCONTACTING  ELECTROOPTIC  DISPLACEMENT  SENSOR  FOR 
PIEZOELECTRIC  ALLY  DRIVEN  ACTIVE  SURFACES* 


S.  A.  Carwin 

Southwest  Research  Institute 
San  Antonio,  Texas 


ABSTRACT 

A  noncontacting  electrooptic  displacement  sensor  for  the  measurement  of  the 
motion  of  a  piezoelectrically  excited  active  wall  is  described.  The  active  wall 
was  constructed  and  studied  as  part  of  a  program  conducted  for  the  Office  of 
Naval  Research  to  investigate  the  drag  reduction  properties  of  an  actively  driven 
surface  in  turbulent  water  flow.  The  sensor  employed  a  two  arm  optical  triangu¬ 
lation  method  with  a  two  element  position  sensitive  detector  to  monitor  the  sur¬ 
face  displacements  from  a  stand-off  distance  of  23  cm  (9  inches) .  Designed  to 
operate  in  either  an  air  or  water  medium,  the  displacement  sensor  demonstrated 
better  than  1.3  pm  (50  p-inch)  resolution  over  displacement  ranges  of  approxi¬ 
mately  1.3  mm  (.05-inch).  The  system  was  implemented  in  two  different  config¬ 
urations:  one  using  an  infrared  LED  light  source  and  the  other  a  visible  red 
helium-neon  laser  source.  In  both  cases,  the  light  source  was  modulated,  and  a 
phase  locked  detector  was  used  to  reject  unwanted  ambient  light.  Spot  size  was 
an  adjustable  parameter  and  was  set  to  1  mm  (.04-inch)  for  this  application.  The 
displacement  monitor  was  articulated  in  three  axes  on  precision  translation 
stages  to  cover  a  5  cm  x  13  cm  (2-in.  x  5-in.)  area.  The  system  was  successfully 
used  to  monitor  accurately  displacements  on  the  order  of  25  is®  (.001-inch)  peak- 
to-peak  on  an  active  wall  driven  by  piezoelectric  transducers  over  an  acoustic 
frequency  range  from  DC  to  150  Hz. 

INTRODUCTION 

This  pape ■  describes  a  noncontacting  electrooptieal  displacement  sensor 
developed  to  mitor  the  motions  of  a  piezoelectrically  excited  active  wall  with 
peak-to-peak  v  ements  on  the  order  of  25  pm.  The  active  wall  was  devised  in 

a  program  conduc.  r  the  Office  of  Naval  Research1,  to  investigate  the  drag 
reduction  prope:tK  *n  actively  driven  surface  in  turbulent  water  flow.  The 

active  portion  of  the  .  was  constructed  by  stretching  a  thin,  diffusely  re¬ 
flecting  mylar  membrane  over  an  array  of  piezoelectric  pushers  spaced  at  an 
interval  of  1.25  mm.  A  cross  sectional  view  of  the  wall  construction  is  given  in 
Figure  1.  Because  each  of  the  peizoelectric  elements  could  be  driven  indepen¬ 
dently,  proper  adjustment  of  the  phase  and  amplitude  of  individual  drive  wave- 


♦Supporteu  by  the  Office  of  Naval  Research  Contract  N00014-82-C-0199 
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the  present  research  was  completed.  Carpenter  and  Garrad^-® 
shown  that  the  Landahl-Kaplan  membrane  model  is  quantita- 
y  incorrect;  however,  they  have  not  computed  the  correct 
■al  stability  curves.  Included  in  Table  II  are  numerical 
.ts  for  both  parallel  and  non-parallel  flow  calculations, 
ion-parallel  calculations^  are  in  better  agreement  with 
riment . 

For  selection  of  dimensional  wave  parameters,  some  physical 
property  must  be  identified.  Since  water  boundary  layers 
/ery  thin,  thickness  is  probably  the  limiting  factor  in 
iary  layer  measurements.  A  reasonable  boundary  layer  thick- 
with  sufficient  resolution  appears  to  be  approximately  2  mm. 
ance  from  the  leading  edge  of  a  flat  plate  for  a  laminar  boun- 
layer  in  a  zero  pressure  gradient  is  related  to  boundary 

r  thickness^  by 

fi/x  =  6.017  (Rex)-1*.  (3) 

nsional  properties  for  T-S  waves  and  a  membrane  are  summarized 
able  III  in  comparison  to  the  boundary-layer  stability  experi- 
s  of  Strazisar  et  al.,2^  for  water. 

On  the  basis  of  information  in  Table  III,  frequency,  wavelength 
wave  speed  can  be  selected  for  the  active  wall;  however,  ampli- 
cannot  be  determined  from  stability  theory.  Normal  velocity 
he  wall  should  be  large  enough  to  influence  the  flow  favorably, 
the  amplitude  should  not  cause  transition.  From  Eq.  (1)  normal 


Lements  which  are  driven  electronically,  in  contrast  to  Kendall's7 
:tive  wall  which  was  driven  mechanically  by  a  cam  and  electric 
5tor.  For  pure  sinusoidal  motion,  coordinates  of  the  wall  are 
Lven  by 

y  =  a  sin  [  a  (x  -  ct)]  (1) 

mere  a  is  the  wave  number,  2tt/\.  In  the  present  application, 
arameters  of  the  wall  motion  are  selected  on  the  basis  of  boun- 
ary  layer  theory. 

Currently,  no  active-wall  boundary-layer  calculations  are 
vailable  for  the  selection  of  parameters,  a,  X  ,  and  c.  However, 
cMurry  et  al.,^-3  in  a  direct  numerical  simulation  have  shown  that 
oth  linear  and  nonlinear  T-S  waves  are  dampened  by  a  vertical  velo- 
ity  component  at  the  wall  which  is  out  phase  by  90°.  With  no 
ther  supporting  data  available,  wavelength  and  wave  speed  were 
elected  from  calculations  of  the  Orr-Sommerf eld  equation. 14 , 15 , 16 
Several  criteria  could  be  used  in  selection  of  the  wave  para- 
eters.  One  criterion  is  the  critical  Reynolds  number,  and  another 
s  derived  from  the  wave  characteristics  of  the  Landahl-Kaplan^7 
embrane  model  that  most  significantly  affect  the  neutral  stability 
urve.  Nondimensional  characteristics  of  the  membrane  model  are 
isted  in  Table  II.  Circular  frequency  w  or  2Ttf,  wave  number  a, 
nd  wave  speed  c  are  nondimensionalized  by  the  freestream  velocity 
oo  and  the  boundary  layer  thickness  6  (thickness  at  99.9%  U^,). 
ave  speed  is  related  to  frequency  and  wavelength  by 

c  =  Xf  =  oj/a.  (2) 
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eir  device  may  only  be  applicable  for  air  flows.  A  summary 
the  wave  characteristics  of  these  active  walls  is  presented 
Table  I. 

Other  devices  have  also  been  successful  in  cancelling  Tollmien- 
ihlichting  waves  and  increasing  transition  Reynolds  number. 

»r  example,  Liepman  and  Nosenchuck10  were  able  to  increase  tran- 
tion  Reynolds  number  in  a  water  tunnel  through  active  control 
:  wall  heating.  Milling^  generated  and  cancelled  T-S  waves 
i  a  water  channel  with  a  pair  of  vibrating  ribbons,  and  Thomas-1-^ 
;rformed  similar  experiments  in  a  wind  tunnel. 

Objectives  of  the  present  active-wall  research  are  much  broader, 
i  addition  to  the  possibility  of  drag  reduction,  the  active  wall 
111  provide  information  on  the  interaction  of  small  traveling 
arface  waves  with  the  boundary  layer.  Such  information  may  identify 
tie  required  wave  characteristics  of  a  passive  compliant  coating 
ar  drag  reduction.  The  results  to  be  described  in  this  paper 
nclude  development  of  performance  criteria,  details  of  the  mechan¬ 
ical  and  electrical  designs,  and  measurement  of  active-wall  perfor- 
ance  characteristics. 

BOUNDARY  LAYER  THEORY 

Conceptually  illustrated  in  Fig.  1,  the  active  wall  consists 
E  a  set  of  actuators  which  are  covered  with  a  flexible  surface, 
hese  actuators  are  driven  sequentially  so  that  a  traveling  surface 
ave  with  amplitude  a,  wavelength  X,  and  phase  speed  c,  is  gener- 
ted.  In  the  present  active  wall,  the  actuators  are  piezoelectric 


INTRODUCTION 


Compliant  coatings  have  been  investigated  over  the  last  25 
years  as  a  means  of  drag  reduction  in  water.  The  first  compliant 
coating  was  developed  by  Kramer, 1' 3  who  speculated  that  dolphin's 
speed  is  enhanced  by  its  pliable  skin.  Since  Kramer  presented 
his  experimental  data  on  the  drag  reducing  properties  of  compliant 
surfaces,  many  theoretical  and  experimental  results  have  been 
published  on  the  physical  mechanisms  involved.  To  date,  these 
results  have  been  rather  inconclusive  according  to  reviews  by 
Bushnell  et  al.3'4  In  many  experiments  drag  increased  with  compliant 
walls,  and  in  those  cases  where  drag  reduction  occurred,  results 
have  not  been  repeatable  by  others. 

A  compliant  coating  of  sufficient  pliability  will  interact 
with  the  flow.  Speculation  is  that  the  boundary  layer  flow  may 
interact  with  a  Kramer-like  coating  so  that  drag  is  reduced. 

However,  very  pliable  coatings  will  form  relatively  large  waves5'6 
that  increase  drag.  Since  the  physics  which  may  allow  drag  to 
be  reduced  is  not  known,  an  active-wall  device  was  designed  and 
constructed  for  a  controlled  investigation  of  the  possibility 
of  drag  reduction  for  small  traveling  surface  waves.  Similar 
research  was  done  by  Kendall, ?  who  used  an  active  wall  to  provide 
insight  into  the  air-sea  interaction  problem  for  a  turbulent  boundary 
layer.  Wehrmann8  successfully  cancelled  Tollmien-Schlichting 
(T-S)  waves  in  an  air  boundary  layer  with  an  active  wall.  An 
interesting  active  wall  concept  based  on  electrostatics  was  developed 
by  Weinstein  and  Balasubramanian^  for  short  wavelengths;  however, 
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fin  all  solid-state  electromechanical  device  has  been  developed 
which  produces  short-wavelength  traveling  surface  waves.  The 
active  wall  was  nominally  designed  for  the  cancellation  of  Tollmien- 
Schlichting  waves  in  a  laminar  boundary  layer  in  water.  The  primary 
purpose  of  the  device  was  determination  of  traveling  surface  wave 
characteristics  which  are  important  for  drag  reduction  with  compliant 
coatings.  Wave  form,  wave  speed,  frequency,  and  amplitude  are 
independently  selectable.  Sinusoidal  wall  motion,  measured  in 
air  with  an  electro-optical  displacement  sensor,  ranged  from  0.6 
to  13  um  over  a  frequency  span  of  10  to  150  Hz.  Roll-off  frequency 
at  the  -3  dB  point  was  120  Hz,  while  the  maximum  phase  speed  was 
166  cm/s.  The  electronic  subsystem  is  a  hybrid  analog/digital 
device,  and  the  active  wall  employs  piezoelectric  ceramic  elements 
to  produce  surface  motion. 

PACS  numbers:  47.80.+V,  47.35.+i,  47.15.Fe,  47.25.Fj 
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an  adjustable  parameter  and  la  controlled  by  the  Q  of  the  5  KHz  bandpass  ampli¬ 
fier,  design  of  the  output  low  pass  filter and  (somewhat)  on  the  selection  of 
carrier  frequency. 

ACTIVE  MALL  MEASUREMENTS 

Actual  measurements  made  on  the  prototype  active  wall  are  given  in  Figures 
6,  7  and  8.  In  Figure  6,  the  amplitude  response  of  an  individual  element  of  the 
active  wall  is  plotted  as  a  function  of  applied  driving  voltage  at  a  fixed  fre¬ 
quency  of  40  Hz .  Figure  7  shows  the  frequency  response  of  the  same  element 
plotted  as  a  function  of  drive  frequency  for  an  applied  voltage  of  150V  peak- 
to— peak.  For  both  Figure  6  and  7  the  raw  data  obtained  from  the  displacement 
monitor  was  corrected  by  the  amplitude  and  frequency  response  calibrations  of 
Figures  4  and  5.  Figure  8  shows  the  longitudinal  phase  response  of  the  active 
wall  when  set  up  to  produce  a  traveling  wave.  The  data  were  taken  by  translating 
the  displacement  monitor  along  the  length  of  the  active  wall  (transversing  from 
element  to  element)  at  a  fixed  standoff  distance,  the  wall  demonstrated  excel¬ 
lent  phase  linearity  over  a  distance  of  approximately  30  mm  (1.2  in.). 

CONCLUSIONS 

In  sununary,  a  methodology  for  producing  noncontacting  displacement  measure- 
ments  through  air  and  water  path  distances  to  surfaces  possessing  very  small 
absolute  displacements  has  been  described.  Actual  measurements  have  been  per¬ 
formed  on  a  piezoelectrical ly  driven  active  well  with  successful  results. 
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Figure  4  is  a  typical  response  plot  of  the  displacement  sensor  taken  with  the 
calibration  fixture*  the  fixture  was  immersed  in  a  transparent  water  tank  for 
this  test,  with  the  reference  surface  set  15  cm  from  the  entry  point  of  the  light 
beam,  the  infrared  LED  was  used  for  the  light  source  and  the  spot  projection  was 
directed  into  the  tank  normal  to  the  top  wall  surface,  the  received  light  was 
collected  through  the  adjacent  tank  side  wall  by  the  receiving  lens  which  was 
oriented  at  45*  to  the  tank  side  wall,  this  rasulted  in  a  watar  path  angle  of 
32*  by  virtue  of  refraction  at  the  water/air  interfaea.  Linearity  and  resolution 
over  the  25  t*a  peak-to-peak  range  is  extrema ly  good,  as  evidenced  in  Figure  4. 
Responsivity  of  the  displacement  sensor  in  this  configuration  was  measured  to  be 
61.58  pm/volt. 

It 


6 

(|im|  * 

-5 

•10 

-0.2  4.1  0  0.1  0.2 

OUTPUT  VOLTAGE 

Figure  4.  Displacement  Calibration  For  Electro-Optical 
Senaor 

Calibration  data  were  taken  in  a  total  of  four  different  configurational  in 
both  air  and  water,  and  with  both  LED  and  laser  light  sources.  The  laser  source 
has  three  distinct  advantages  and  two  distinct  disadvantages  for  this  applica¬ 
tion.  First,  the  laser  provides  a  much  more  intense  source  of  light  and  operates 
at  a  wavelength  which  has  little  loss  in  water.  The  LED  source,  at  933  m,  is 
attenuated  by  the  absorption  band  in  water  near  944  nm,  presenting  a  signal 
strength  problem  for  long  water  path  distances.  Second,  the  laser  can  be  focused 
to  a  much  smaller  spot  on  the  target  surface,  providing  superior  spatial  resolu¬ 
tion.  The  LED  source  cannot  be  focused  efficiently  to  a  spot  much  smaller  than 
the  physical  diameter  of  the  active  area  without  incurring  significant  losses  in 
intensity.  Finally,  the  He-Ne  laser  operates  in  the  visible  red  (as  opposed  to 
IR  for  the  LED),  making  alignment  of  the  system  an  easier  task. 

The  disadvantages  of  using  a  laser  are  expense  and  stability.  Most  commer- 
cially  available  lasers  in  this  power  range,  although  relatively  inexpensive  in 
themselves,  cannot  be  directly  modulated.  An  aoousto— optical  modulator  was  used 
to  provide  the  required  5  kHz  subcarrier,  and  the  cost  of  the  modulator  was  more 
than  that  of  the  laser  itself.  Secondly,  the  laser  must  have  a  stabilized  output 
relatively  free  from  the  noise  components  associated  with  the  Intermodulation 
products  of  multiple  optical  frequencies.  These  products  commonly  occur  in  the 
frequency  range  of  1  to  100  kHz,  which  is  in-band  to  the  5  kHz  subcarrier  fre¬ 
quency.  The  noise  components  appear  as  sidebands  about  the  5  kHz  subcarrier 
frequency  and  thus  cannot  be  filtered  easily.  The  analog  division  process  re¬ 
moves  most  of  the  amplitude  modulated  noise  components,  but  improvements  in  sig¬ 
nal  to  noise  ratios  can  be  gained  through  the  use  of  stabilised  lasers. 

The  frequency  response  of  the  displacement  monitor  was  measured  by  inserting 
an  analog  multiplier  between  the  detector  preamplifier  and  the  remainder  of  the 
signal  processing  electronics.  Light  input  to  the  detectors  was  held  constant 
while  the  frequency  of  an  external  oscillator  connected  to  the  control  input  of 
the  multiplier  was  varied.  The  frequency  response  curve  is  given  in  Figure  5  and 
shows  reasonably  flat  response  out  to  about  60  Hz  with  a  roll  off  of  approxi¬ 
mately  8  dB/ octave  after  the  breakpoint.  Response  of  the  displacement  sensor  is 


The  sum  of  ths  two  detector  half  outputs  is  a  constant  with  raspset  to  spot 
position  and  contains  spot  intensity  and  diameter  information.  The  difference 
function  is  also  influenced  directly  by  spot  intensity  and  diameter  wfiich  is  un¬ 
desirable,  as  errors  in  the  displacement  measurement  will  occur  if  factors  in¬ 
fluencing  spot  intensity  or  detector  sensitivity  (o.g.  LED  or  laser  brightness, 
target  surface  reflectivity,  water  path  attenuation,  temperature,  etc.)  are  vari¬ 
able.  As  the  intensity  and  diameter  data  appear  in  both  the  difference  and  sum 
function,  dividing  the  difference  by  the  sum  cancels  errors  associated  with  vari¬ 
able  intensity  or  sensitivity.  This  method  assumes  uniform  intensity  within  the 
spot,  and  no  difficulties  were  encountered  with  nonuniform  illumination  or  re¬ 
flectance  in  this  application.  For  details  on  a  distance  measurement  system 
based  on  a  similar  principle  but  with  provisions  to  cancel  errors  associated  with 
nonuniformly  reflective  surfaces,  the  interested  reader  is  referred  to  the  system 
developed  for  profiling  roadway  surfaces  [2] . 

A  block  diagram  of  the  signal  processing  electronics  for  the  displacement 
sensor  is  given  in  Figure  3.  Design  of  the  circuitry  is  straightforward  and 
models  the  foregoing  signal  processing  algorithm.  A  S  kHz  crystal  controlled 
oscillator  is  used  to  drive  the  light  modulator  and  to  act  as  reference  for  the 
phase  sensitive  detectors.  The  outputs  of  the  two  detector  halves  are  individ¬ 
ually  preamplified  and  filtered  through  active  5  kHz  bandpass  filters.  DC  re¬ 
storation  is  accomplished  through  phase  sensitive  detectors  and  subsequent  ampli¬ 
fication  by  DC  amplifiers.  The  DC  voltages  thus  obtained  (denoted  by  El  and  E2) 
are  proportional  to  the  light  values  impinging  on  each  detector  half.  El  and  E2 
are  then  processed  by  the  differential  and  summation  amplifiers  and  the  resultant 
difference  and  sum  are  fed  into  the  Y  and  X  inputs  of  the  analog  divider. 

The  quotient  is  (E1-E1 )/(E1+E2)  and  is  taken  as  the  displacement  signal  after 
high  frequency  noise  component  are  removed  by  an  active  low-pass  filter.  The 
signal  processing  circuitry  was  packaged  in  a  compact  (12  cm  x  18  cm  x  25  cm) 
enclosure  and  was  completely  self-contained,  Bie  light  source,  lens  systems,  and 
detector  assembly  were  mounted  on  a  precision  XYZ  fixture  for  positioning  over 
selected  areas  of  the  active  wall. 

Specification  for  some  of  the  component  parts  implemented  in  the  laboratory 
system  are  as  follows t 

LED  Source:  Type  TIES 12  (Texas  Instrument) 

0.91  nn  (0.036-inch)  diameter 
50  ow  (  933  nm 

LED  Projection  lens  72  mm  focal  length 

Assembly:  SO  ns  diameter 

Detector:  Type  NN  Spot-2D  (United  Detector 

Technology)  two  element  discreet 
active  surface:  2.54  x  2.67  mm  (0.1  x  0.105-in.) 

Detector  Lens  47  ns  focal  length 

Assembly :  45  am  diameter 

Laser  source:  Helium-Neon  Type 

Spectra- Physics 

model  138,  1.0  mw  9  633  nm 

Laser  Modulator:  Acoustooptic  Type 

Anderson  labs  Model  PLM-5VS 
40  mHz  center  acoustic  frequency 
100%  square  wave  modulation  0  5  kHz 

SENSOR  CALIBRATION  AND  FREQUENCY  RESPONSE 

The  electrooptic  displacement  sensor  was  calibrated  with  a  fixture  con¬ 
structed  specifically  for  the  task.  The  calibration  fixture  was  fabricated  by 
using  a  precision  micrometer  and  ball  configuration  in  conjunction  with  a  10:1 
mechanical  reduction  arm  to  provide  precise  control  over  the  reference  surface. 
Thus  0.010-inch  (254  um)  of  travel  on  the  micrometer  head  (which  was  readable  to 
0.0001-inch  or  2.54  um)  produced  0.001-inch  (25.4  um)  of  travel  at  the  reference 
surface.  Precision  ball  bearings  were  used  for  the  fulcrum. 
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Figure  2.  Optical  Displacement  Measuring  Technique 


Figure  3.  Electro-Optical  Displacement  Sensor 
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forms  provided  stationary  or  traveling  surface  waves  programmable  in  direction, 
amplitude,  frequency,  and  in  the  case  of  traveling  waves,  wave  speed.  Details  of 
the  entire  active  wall  program  are  being  presented  elsewhere  [1]  and  the  remaind¬ 
er  of  this  paper  describes  the  electrooptieal  displacement  sensor  and  subsequent 
measurements  performed  on  the  active  wall. 


Figure  1.  Illustration  of  Piezoelectrleally  Driven 
Active  Wall. 


PRINCIPLES  OF  OPERATION 

The  methodology  used  for  the  measurement  technique  is  an  adaptation  of  a 
technique  developed  by  the  author  to  profile  roadway  surfaces  from  a  moving 
vehicle  (2]  and  is  illustrated  in  Figure  2.  Two  variations  of  the  method  can  be 
implemented  by  choice  of  light  source:  a  1  mw  Helium-Neon  laser  operating  in  the 
visible  red  at  633  nm  or  a  high  power  (100  mw)  light  emitting  diode  (LED)  opera¬ 
ting  in  the  near  infrared  at  933  nm.  In  both  cases  the  light  source  is  100% 
amplitude  modulated  at  a  frequency  of  S  kHz  for  use  with  phase  sensitive  detec¬ 
tion  so  that  a  high  signal-to-noise  ratio  and  rejection  of  ambient  light  can  be 
achieved.  For  the  infrared  version,  the  LED  is  modulated  directly  by  the  drive 
current,  and  for  the  visible  red  version,  an  acousto-optic  modulator  is  used  to 
modulate  the  laser.  The  light  source  is  projected  onto  the  target  surface  by  a 
lens  assembly  oriented  normal  to  the  target  surface  to  form  a  spot  lnmt  in  diam¬ 
eter,  a  spot  size  small  enough  to  resolve  individual  piezoelectric  elements. 

A  portion  of  the  light  scattered  from  the  diffusely  reflecting  target  sur¬ 
face  is  collected  by  a  receiving  lens  assembly  oriented  at  45*  to  the  target  sur¬ 
face.  The  lens  assembly  images  the  illuminated  spot  onto  the  center  of  a  two 
element  PIN  photodiode  detector  with  unity  magnification.  This  target  surface 
position,  denoted  as  POSITION  1  in  Figure  2,  is  the  reference  (or  zero),  position 
about  which  displacements  are  measured.  At  this  position,  the  spot  image  is 
exactly  centered  on  the  two  detector  halves  thereby  causing  the  outputs  of  the 
two  halves  to  be  equal.  By  virtue  of  the  45*  geometry,  a  positive  or  negative 
displacement  of  the  target  surface  from  the  reference  position  causes  a  corres¬ 
ponding  lateral  shift  of  the  spot  image  on  the  face  of  the  dual  detector.  An 
example  of  a  negative  displacement  is  shown  as  POSITION  2  and  is  represented  by 
the  dashed  lines  for  the  target  surface,  light  path,  and  spot  image  on  the  dual 
detector.  Shifting  the  image  of  the  spot  on  the  face  of  the  detector  causes  more 
area  of  the  spot  to  fall  on  one  half  of  the  detector  than  the  other,  thus  pro¬ 
ducing  unequal  outputs.  The  displacement  signal  is  extracted  by  computing  the 
difference  and  sum  of  the  outputs  of  the  two  detector  halves,  then  dividing  the 
difference  by  the  sum.  The  difference  between  the  two  detector  halves  as  the 
image  of  the  spot  is  translated  from  one  detector  half  to  the  other,  (beginning 
with  the  spot  completely  one  half) ,  is  simply  the  difference  in  areas  of  a  circle 
divided  by  a  chord.  The  function  is  *S*  shaped,  with  extremely  good  linearity 
for  small  displacements  about  center. 


v  =  dy/dt  =  -aw  cos  [ct(x  -  ct)). 


(4) 


Then,  amplitude  is 


a/6  =VT  (av/Uo  )  /  (w6  / U«  ) 


(5) 


where  ov  is  the  standard  deviation  of  the  normal  velocity  at  the 
wall.  The  normal  relative  intensity,  ov/Uoo,  at  the  wall  should 
be  higher  than  the  relative  turbulence  intensity  which  is  0.1% 
in  the  better  water  tunnels . 20 • 21 • For  an  adequate  rms  signal- 
to-noise  ratio,  ov/Uoo  should  be  an  order  of  magnitude  larger; 
consequently,  from  an  instrumentation  viewpoint  av/V00  would  be 
measured  in  the  flow  with  an  accuracy  of  10%  for  an  rms  signal- 
to-noise  ratio  of  10  .  With  a  normal  velocity  component  of  1% 
and  a  nondimensional  frequency  of  0.5,  nondimensional  amplitude 
should  be  0.028.  For  a  boundary  layer  thickness  of  2  mm,  amplitude 
should  then  be  57  urn.  However,  on  the  basis  of  other  design  con¬ 
siderations  which  are  to  be  described  subsequently,  a  design  goal 
of  25  pm  was  selected. 

Some  experimental  evidence  is  available  on  the  required 
amplitude  of  the  normal  wall  motion.  Wehrmann®  was  able  to  cancel 
T-S  waves  in  a  wind  tunnel  with  an  rms  normal  velocity  of  3  X  10“3% 
which  is  an  order  of  magnitude  smaller  than  the  relative  turbulence 
intensity  of  a  good  wind  tunnel. 23  in  the  strip  heater  experiments 
of  Liepmann  et  al.22,  an  equivalence  between  wall  heating  and 
wall  motion  was  derived  for  a  water  boundary  layer.  For  their 
experiments,  the  rms  normal  wall  motion  was  related  to  the  maximum 
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wall  overheat,  AT,  in  °C  by 

ov/Uoo  =  10"4  | A  T  | .  (6) 

The  overheat  was  typically  3°c;  therefore,  the  equivalent  rms 
wall  motion  was  0.03%  which  was  the  same  order  of  magnitude  as 
the  0.04%  longitudinal  turbulence  intensity  of  the  water  tunnel. 22 

Eight  driving  elements  per  wavelength  were  arbitrarily  selected 
as  adequate  for  generation  of  a  sine  wave.  WehrmannS  used  ten 
in  his  active  wall  while  Kendall?  claimed  little  distortion  with 
four.  On  the  basis  of  cost  and  the  rationale  of  the  previous 
discussion,  the  design  criteria  summarized  in  Table  IV  were  chosen. 
Comparison  of  Tables  III  and  IV  indicates  that  the  design  criteria 
meet  or  exceed  conditions  likely  to  be  encountered  in  a  water 
flow  experiment. 

II.  PRINCIPLES  OF  OPERATION 
A.  Piezoelectric  Actuators 

One  of  the  major  design  goals  for  an  active  wall  is  that 
amplitude,  frequency,  and  speed  of  the  surface  wave  should  be 
independently  controllable.  Most  schemes  for  generating  active- 
wall  motion  permit  arbitrary  control  of  only  one  or  two  of  these 
parameters.  For  example,  only  wave  speed  could  be  varied  in 
Kendall's?  experiments 

For  the  present  active  wall,  piezoelectric  ceramic  materials 
were  selected  as  the  most  suitable  electromechanical  conversion 
actuators.  A  schematic  cross  section  of  a  piezoelectrically 
actuated  active  wall  is  shown  in  Fig.  1  in  which  a  number  of  piezo¬ 
electric  elements  with  long  dimensions  vertical  are  standing  on  a 


rigid  base.  Each  element  has  a  different  voltage  applied  to  it 
so  that  different  degrees  of  elongation  (greatly  exaggerated  in 
the  figure)  exist  at  each  element.  A  thin  membrane  attached  to 
the  tops  of  all  of  the  piezoelectric  elements  provides  for  inter¬ 
polation  to  produce  relatively  smooth  curvature  of  the  wall  surface 
and,  of  course,  prevents  fluid  penetration  into  the  volume  occupied 
by  the  piezoelectric  elements. 

In  operation,  all  of  the  piezoelectric  elements  are  driven 
by  electrical  signals  of  the  same  amplitude  and  frequency;  however, 
the  signal  applied  to  each  successive  element  proceeding  from 
left  to  right  in  Fig.  1  is  delayed  in  phase  by  a  progressively' 
increasing  amount  so  that  a  sinusoidal  disturbance  is  produced 
on  the  active-wall  surface.  Additionally,  at  a  selected  time 
interval,  the  instantaneous  signal  value  on  each  element  is  shifted 
to  the  next  element  to  the  right  so  that  the  sinusoidal  disturbance 
propagates  from  left  to  right. 

A  piezoelectric  material  changes  physical  dimensions  in  response 
to  an  applied  electric  field;  conversely,  when  a  piezoelectric 
material  is  mechanically  stressed,  it  generates  a  proportionate 
electrical  signal.  Hence,  piezoelectric  materials  may  be  used 
as  electrical-to-mechanical  or  mechanical-to-electrical  transducers. 
For  the  active  wall,  the  electrical-to-mechanical  conversion  capa¬ 
bility  is  utilized.  Figure  2  illustrates  a  rectangular  piezo¬ 
electric  ceramic  slab  having  dimensional  proportions  that  dictate 
operation  in  the  longitudinal-extension  mode;  that  is,  length  is 
at  least  three  times  width  ,  and  width  is  substantially  greater 


than  thickness.  Further,  the  two  large  faces  (defined  by  length 
and  width)  are  covered  with  metallic  electrodes  between  which  an 
electrical  drive  voltage  is  applied. 

In  Fig.  2  a  system  of  coordinate  axes  is  defined,  and  these 
axes  are  labeled.  The  direction  of  polarization  is  designated 
axis  3,  and  the  principal  mechanical  axis  is  designated  1.  The 
remaining  axis  in  the  mutually  perpendicular  triplet  is  designated  2. 
As  indicated  in  Fig.  2,  when  an  electrical  drive  voltage  is  applied 
between  electroded  faces  of  the  piezoelectric  element,  elongation 
(or  contraction,  which  is  negative  elongation)  occurs  along  the 
principal  mechanical  axis.  Whether  elongation  is  positive  or 
negative  depends  upon  polarity  of  the  applied  drive  voltage. 

Relationships  between  electrical  drive  and  mechanical  response 
are  described  by  a  series  of  piezoelectric  constants  characteristic 
of  the  particular  material  employed.  In  the  case  of  a  longitudinal- 
extension  element  such  as  illustrated  in  Fig.  2,  the  constant 
of  proportionality  is  the  piezoelectric  strain  constant  d3i  which 
is  the  ratio  of  strain  (i.e.,  ratio  of  elongation  to  length)  to 
impressed  electric  field  (i.e.,  ratio  of  applied  voltage  to  thick¬ 
ness).  The  impressed  electric  field  is  in  the  3  direction  and 
mechanical  response  is  in  the  1  direction;  hence,  the  subscript  31 . 
Expressed  mathematically,  the  relationship^ , 25  between  mechanical 
response  and  impressed  electric  field  is 

o  =  d3^e  (7) 

where  0  is  strain,  £  is  electric  field,  and  d3i  is  piezoelectric 
strain  constant. 
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By  definition, 


a  =  A£/£  (8) 

and 

E  =  V/t  (9) 

where  A£  is  elongation,  £  is  element  length,  V  is  applied  voltage, 
and  t  is  element  thickness.  Prom  Equations  (7),  (8),  and  (9) 

A£  =  (d31£V)/t  .  (10) 

From  Equation  (10)  piezoelectric  strain  constant,  element  length, 
and  applied  voltage  should  be  maximized  while  element  thickness 
should  be  minimized  for  the  required  elongation,  but  practical 
constraints  on  piezoelectric  element  dimensions  and  on  magnitude 
of  applied  voltage  exist.  Further,  only  a  limited  range  of  piezo¬ 
electric  strain  constants  is  available. 

Current  technology  limits  the  size  of  sheet  piezoelectric 
materials  to  a  maximum  length  of  76  mm  and  a  minimum  thickness 
of  0.25  mm.  Thus,  the  limiting  values  of  length  and  thickness 
(76  mm  and  0.25  mm,  respectively)  and  a  width  of  25  mm  (i.e.  one- 
third  of  the  length  to  insure  longitudinal-extension  mode  of  oper¬ 
ation)  were  selected  for  element  dimensions.  Further,  a  lead 
zirconate  titanate  (PZT)  ceramic,  PZT-5H  from  Vernitron,  Inc., 
Bedford,  Ohio,  was  chosed  because  it  has  the  highest  available 
piezoelectric  strain  constant  (d3i)  with  a  value  of  274  pm/V. 

With  these  parameters  and  a  peak  applied  voltage  of  300  V,  peak 
elongation  or  amplitude  was  estimated  to  be  25  um  from  Eq.  (10). 


B.  Mechanical  Assembly 

The  active-wall  model  was  built  in  cells  which  were  fabricated 
as  shown  in  Fig.  3.  Each  cell  contained  two  piezoelectric  elements 
with  dimensions  76-mm  long  X  25-mm  wide  X  0.27-mm  thick.  Two 
such  elements  were  placed  side  by  side  for  an  assembly  that  was 
51-mm  wide.  A  0.03-mm  thick  metal-foil  electrode  with  an  integral 
electrical  connection  tab  was  placed  on  each  side  of  the  pair 
of  piezoelectric  ceramic  plates  with  the  connecting  tabs  from 
opposite  faces  extending  from  opposite  sides  of  the  assembly. 

A  silicone  rubber  insulator  0.46-mm  thick  was  placed  on  each  side 
of  the  assembly  to  complete  an  individual  cell.  The  silicone 
rubber  insulators  between  piezoelectric  elements  provided  not 
only  electrical  insulation,  but  also  a  compliant  interface  to 
insure  independence  of  operation  of  adjacent  piezoelectric  elements. 
Total  cell  thickness  was  1. 25-mm.  Elongation  is  parallel  to  the 
75-mm  dimension  of  the  cell,  and  fluid  flow  is  parallel  to  the 
thickness  of  the  cell  as  shown  in  Fig.  3. 

Individual  cells  were  stacked  with  large  face  to  large  face 
in  a  box-like  container  as  shown  in  Fig.  4.  Because  each  cell 
had  a  layer  of  silicone  rubber  on  each  side,  in  the  assembly  two 
thicknesses  of  silicone  rubber  were  located  between  pairs  of  piezo¬ 
electric  elements.  Directions  of  elongation  and  fluid  flow  relative 
to  the  overall  structure  are  also  shown  in  Fig.  4.  The  box  was 
constructed  of  Plexiglas. 

The  top  surface,  which  acted  as  the  interface  between  the 
flowing  fluid  and  the  cavity  containing  the  piezoelectric  elements, 


was  a  thin  film  of  adhesive-backed,  shrinkable  plastic  film  (model 
airplane  covering)  that  had  an  opaque  white  finish.  The  film 
was  adhered  to  the  top  edges  of  the  plastic  box  and  then  was  ther¬ 
mally  shrunk  to  make  it  contact  the  tops  of  the  elements.  The 
resulting  active-wall  area  was  51  mm  wide  (transverse  to  the  fluid- 
flow  direction)  and  40  mm  long  in  the  fluid-flow  direction. 

C.  Electronic  Control 

1.  Master  Control  Circuit 

The  basic  system  for  generating  drive  signals  for  the  individual 
piezoelectric  actuator  cells  involves  phase  shift  of  a  sine  wave 
while  its  other  characteristics  are  maintained  constant.  Phase 
shift  was  achieved  with  digital  shift  registers.  The  low-level 
signal-generating  circuits  were  divided  into  two  functional  cate¬ 
gories:  (1)  the  master  control  circuit  and  (2)  the  individual 

channel  circuit.  The  master  control  circuit  provided  an  interface 
between  the  analog  output  of  an  external  laboratory  signal  generator 
and  a  digital  signal  that  could  be  passed  from  one  channel  to 
the  next.  The  control  system  contained  32  identical  channel  cir¬ 
cuits  in  cascade. 

Figure  5  illustrates  the  master  control  circuit  in  block 
diagram  form.  The  input  signal  is  derived  from  a  commercially 
available  laboratory  function  generator  that  provides  sinusoidal 
output  with  continuously  variable  amplitude  and  frequency.  The 
input  signal  is  converted  periodically  to  an  8-bit  offset-binary 
digital  word  whenever  an  electronic  command  is  received  from  the 
timing  generator.  The  analog-to-digital  converter  block  also 


includes  a  sample-and-hold  circuit  so  that  the  input  signal  does 
not  change  during  the  conversion  period. 

System  timing  was  controlled  by  a  variable-frequency  square- 
wave  oscillator  that  was  an  integral  part  of  the  master  control 
circuit.  By  establishing  system  timing,  this  oscillator  controlled 
wave  speed  of  the  traveling  wave  on  the  active-wall  surface. 

That  is,  a  given  sample  of  the  input  signal  advanced  from  one 
piezoelectric  element  in  the  active-wall  array  to  the  next  every 
time  a  clock  pulse  occurred.  Wave  speed,  then,  was  the  distance 
that  the  wave  advanced  along  the  active-wall  surface  (i.e.,  the 
cell-to-cell  spacing  a*  of  1.25  mm)  divided  by  the  time  interval 
Ars  between  clock  pulses.  Since  the  highest  frequency  that  the 
power  amplifier  could  safely  accommodate  was  150  Hz  and  the  minimum 
number  of  samples  per  cycle  was  8,  the  minimum  design  sampling 
interval  was  833  us,  and  the  maximum  wave  speed  was  then  1.5  m/s. 
Whenever  the  sampling  interval  was  longer  than  the  833-us  minimum, 
wave  speed  was  reduced  proportionately. 

When  the  output  of  the  variable-frequency  oscillator  changed 
logic  state  in  a  selected  direction,  the  timing  generator  issued 
a  convert  command  to  the  analog-to-digital  converter.  At  this 
point  the  input  signal  was  sampled  and  was  converted  to  digital 
form.  The  timing  generator  then  waited  until  analog-to-digital 
conversion  was  complete.  After  a  valid  digital  word  was  established 
on  the  8-line  digital  output,  a  simultaneous  clock  pulse  was  issued 
by  way  of  a  buffer  to  all  32  individual  channel  circuits. 

Frequencies  in  the  range  of  10  Hz  to  150  Hz  could  be  accepted 

by  the  electronic  system.  Going  beyond  the  upper  limit  would  result 
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in  overstress  of  the  transistors  in  the  power  amplifier.  Wave 
speed  could  range  from  10  cm/s  to  a  maximum  that  was  numerically 
equal  to  the  input  frequency.  For  example,  at  an  input  frequency 
of  100  Hz  the  maximum  wave  speed  was  100  cm/s.  This  upper  limit 
insured  that  waves  on  the  active  surface  would  contain  at  least 
8  piezoelectric  cells  per  wavelength.  If  lower  resolution  could 
be  tolerated,  then  the  wave  speed  limit  could  be  exceeded  without 
danger  of  damage  to  the  circuit. 

2.  Individual  Channel  Circuit 

The  circuit  associated  with  an  individual  channel  of  the 
active-wall  driver  is  illustrated  in  the  block  diagram  of  Fig.  6. 

The  8-line  digital  input  originates  at  the  output  of  the  master 
control  circuit  analog-to-digital  converter  (for  channel  1  only) 
or  at  the  output  of  the  previous  channel  (for  channels  2  through 
32).  Each  channel  incorporated  a  single-stage  digital  shift  register 
with  the  capability  of  storing  and  simultaneously  transferring 
8  lines  of  digital  information.  As  previously  noted,  the  output 
of  the  shift  register  in  a  particular  channel  became  the  digital 
input  for  the  subsequent  channel.  Upon  receipt  of  a  clock  pulse, 
which  entered  an  individual  channel  on  the  clock  input  line  and 
propagated  to  the  next  channel  by  way  of  the  clock  output  line, 
the  8-bit  digital  word  present  on  the  8-line  digital  input  was 
stored  in  the  shift  register,  and  the  digital  word  that  was  in 
the  shift  register  prior  to  the  clock  pulse  was  stored  in  the 
next  channel  shift  register  and  so  on.  Thus,  a  given  sample  of 
the  analog  input  signal  moved  without  distortion  or  attenuation 
from  channel  to  channel  when  commanded  to  do  so  by  clock  pulses. 


When  a  digital  word  appeared  at  the  output  of  a  particular 
shift-register  stage,  a  proportionate  analog  current  was  produced 
by  a  digital-to-analog  converter.  This  current  was  converted 
to  a  voltage  by  a  buffer  amplifier,  and  the  resulting  voltage 
was  applied  to  a  power  amplifier  which,  in  turn,  drove  the  piezo¬ 
electric  ceramic  element  in  the  active-wall  cell  corresponding 
to  that  channel.  Thus,  as  sequential  clock  pulses  occurred,  the 
analog  signal  impressed  on  a  selected  piezoelectric  element  replicated 
the  master  control  circuit  input  signal  delayed  by  an  integral 
number  of  clock-pulse  intervals. 

At  frequencies  well  below  primary  resonance,  electrical  charac¬ 
teristics  of  a  piezoelectric  element  are  essentially  equal  to 
those  of  a  capacitor.  Driving  any  kind  of  reactive  load  with 
a  power  amplifier  results  in  a  nonlinear  load  line  for  the  power 
semiconductor  devices  in  the  output  stage.  Consequently,  current 
through  the  power  device  is  high  when  voltage  across  the  device 
is  also  high;  therefore,  substantial  instantaneous  heating  occurs 
which  must  be  dissipated.  Permissible  combinations  of  voltage 
across  a  transistor  and  current  through  it  are  illustrated  graphi¬ 
cally  in  Fig.  7  by  the  safe  operating  areas  (SOA)  from  manufacturers' 
specifications  for  bipolar  transistors  and  vertical  metal-oxide- 
silicon  field-effect  transistors  (VMOSFET)  with  the  same  voltage 
rating.  The  load  line  for  a  0.44-yF  capacitor,  which  simulates 
the  total  capacitance  of  an  active-wall  cell,  driven  at  150  Hz 
by  a  600-V  peak-to-peak  sinusoid,  is  also  plotted  in  Fig.  7. 

As  the  figure  indicates,  the  VMOSFET  is  the  most  suitable  power 

amplifier  in  the  present  application. 

33 


V  V  *j 


*•  .% 


,  •  ■  *  4  .  4  '  .  *  ■  4  4  •  |  «  4  ■  4  •  •  •  *  t  4  < 

»  •  -  ’  •  WJ *.*  V  V-V.V- V_ * .•  V  V  ■; '/.-.V 


To  develop  the  greatest  possible  voltage  across  the  capacitive 
load  from  minimum  power-supply  voltage  requires  use  of  a  bridge 
amplifier  configuration  such  as  that  illustrated  in  Fig.  8  in 
which  details  of  the  bias  circuits  have  been  eliminated  for  clarity 
The  bridge  circuit  was  implemented  with  complementary- symmetry 
transistors  so  that  all  of  the  bias  networks  in  Fig.  8  are  referred 
to  the  power  supplies,  and  the  drive  voltages  required  by  the 
upper  and  lower  transistors  are  relatively  low  in  amplitude  and 
equal.  The  transistors  ultimately  used  in  the  individual  channel 
circuits  of  the  active-wall  driver  had  voltage  ratings  of  350  V 
which  would  have  facilitated  developing  700  V  peak-to-peak  across 
the  capacitive  load.  However,  the  maximum  voltage  available  from 
suitable  commercial  power  supplies  was  300  V,  and  this  placed 
a  600-V  peak-to-peak  limitation  on  the  signal  that  could  be  devel¬ 
oped  across  the  capacitive  load. 

In  the  circuit  shown  in  Fig.  8  a  positive-swinging  voltage 
applied  to  the  gate  (G)  of  p-channel  VMOSFET  Ql  causes  that  tran¬ 
sistor  to  cut  off,  while  the  same  positive-going  voltage  applied 
to  the  gate  of  n-channel  VMOSFET  Q3  causes  that  transistor  to 
turn  on.  A  signal  of  the  same  amplitude,  but  inverted  (i.e., 
shifted  180°  in  phase),  is  applied  simultaneously  to  the  gates 
of  VMOSFET' s  Q2  and  Q4  producing  similar  action  so  that  Ql  and 
Q4  are  simultaneously  turned  off  while  Q2  and  Q3  are  turned  on. 

In  this  case  the  right-hand  end  of  the  capacitive  load  is  pulled 
to  the  positive  supply  by  Q2  while  the  left-hand  end  of  the  capa¬ 
citive  load  is  pulled  to  the  negative  supply  by  Q3 .  During  the 


>osite  half -cycle  of  the  sinusoidal  input  voltage,  roles  of 
;  VMOSFET  pairs  reverse;  that  is,  Ql  and  Q4  are  turned  on  thereby 
Lling  the  left-hand  end  of  the  load  to  the  positive  supply  and 
s  right-hand  end  of  the  load  to  the  negative  supply.  At  the 
tie  time,  VMOSFET' s  Q2  and  Q3  are  turned  off  so  that  no  current 
ows  through  them.  Thus,  voltage  across  the  load  varies  sinusoid- 
ly  from  +2  V  (where  V  is  the  voltage  of  either  the  positive 
negative  power  supply)  to  -2  V  during  the  opposite  half -cycle 
the  input  signal;  consequently,  peak-to-peak  voltage  appearing 
ross  the  capacitive  load  is  4  V.  In  the  individual  channel 
wer  amplifiers  constructed  for  this  program,  the  value  of  V 
s  150  V,  so  that  the  maximum  voltage  developed  across  the  load 
s  600  V  peak-to-peak.  A  low-voltage  zener  diode  was  incorporated 
i  each  gate  circuit  so  that  the  gate-to-source  voltage  could 
>t  exceed  rated  value  during  power-supply  turn  on. 

VMOSFET' s  are  enhancement-mode  devices  which  means  that  when 
le  gate  and  source  are  at  the  same  potential  (i.e.,  zero  bias), 

>  drain  current  flows  through  the  transistor.  As  the  gate-to-source 
>ltage  is  increased  in  the  same  polarity  as  the  drain-to-source 
iltage,  a  threshold  is  reached  at  which  the  transistor  suddenly 
:arts  conducting  and  transconductance  (i.e.,  gain)  increases  drama- 
cally.  The  latter  characteristic  implies  that  small  additional 
icreases  in  gate-to-source  voltage  will  cause  large  changes  in 
•ain  current;  hence,  unless  gate-to-source  voltage  is  carefully 
mtrolled,  destructive  currents  may  flow  through  the  transistor, 
le  problem  is  compounded  because  substantial  device-to-device 


35 


iriation  occurs  in  threshold-voltage  and  transconductance  charac- 
sristics.  Hence,  a  fixed  bias  network  may  leave  one  device  turned 
: f ,  but  may  cause  destructive  currents  to  flow  through  another 
;vice  from  the  same  manufacturing  lot.  These  problems  were  over- 
>me  by  providing  an  adjustment  in  each  bias  circuit  and  developing 
safe  procedure  for  setting  the  operating  point  of  each  transistor 
1  the  bridge  power  amplifier. 

El.  ACTIVE  WALL  MEASUREMENTS 

Motion  of  the  active  wall  in  air  was  measured  with  an  electro- 
ptical  displacement  sensor  which  employed  a  two-arm  optical  triangu- 
ation  method  and  a  two-element  position-sensitive  detector  to 
onitor  displacements  at  a  point  from  a  stand-off  distance  of 
3  cm.  Sensor  response  was  very  linear  over  the  range  of  +  13  pm 
ith  resolution  less  than  0.5  pm  and  sensitivity  of  61.5  pm/V. 
requency  response  had  a  low-pass  characteristic  with  a  -3  dB 
oll-off  frequency  of  90  Hz.  Since  the  active  wall  was  designed 
or  an  upper  frequency  limit  of  150  Hz,  measurements  of  the  active 
all  were  corrected  with  amplitude-response  data.  Additional 
etails  are  discussed  by  Cerwin.^6 

The  active  wall  was  mounted  on  an  x-y  translator  with  micrometer 
djustments  so  that  the  surface  of  the  active  wall  could  be  located 
ormally  and  longitudinally  relative  to  the  displacement  sensor 
ith  a  resolution  of  2.5  pm  over  the  24. 5- mm  travel  ranges  of 
he  micrometers.  The  active  wall  was  calibrated  for  amplitude 
nd  frequency  response  at  a  station  between  the  first  and  second 
ow  of  elements.  At  a  frequency  of  40  Hz,  the  input  voltage  ampli- 
ude  to  the  elements  was  varied  between  10  and  150  V  with  resulting 


acement  amplitude  variation  between  0.62  and  12.9  ym,  respec- 
y.  Distortion  of  the  sine  wave  impressed  on  the  elements 
red  to  occur  between  80  and  90 -V  amplitude.  Calibration 
ts  are  plotted  in  Fig.  9  and  the  amplitude  response  curve 
he  active  wall  is  presented  in  Fig.  10  for  two  input  voltage 
tudes  normalized  at  10  Hz.  The  roll-off  frequency  at  -3  dB 
0  Hz  for  the  50-V  input  amplitude. 

From  data  in  Fig.  9  and  physical  dimensions  of  the  ceramic 
:nts ,  the  piezoelectric  strain  constant  d3i  was  computed  from 
10).  The  mean  value  of  d3i  for  a  single  element  over  a  voltage 
;  of  10  to  140-V  peak  at  40  Hz  was  259  pm/V  at  a  standard 
ition  of  54  pm/V  which  is  within  5.5%  of  the  manufacturer's 
Lfication.  Displacement  measurements  were  also  taken  at  ten 
srent  elements  at  100- V  peak  and  40  Hz.  The  d3i  for  these 
arements  was  264  pm/V  with  a  standard  deviation  of  51  pm/V. 
Since  the  displacement  sensor  is  a  point  measurement  device, 

2  relation  between  longitudinal  locations  on  the  active  wall 
neasured  to  demonstrate  that  the  surface  waves  were  traveling, 
phase  shift  <p  in  degrees  which  is  designed  into  the  active 
is  given  by 


<p  =  36 Of  ( At  S/Ax  )x  (11) 

b  Ax  is  the  spacing  between  elements,  Axs  is  the  sample  interval 
le  clock  in  the  master  control  circuit  of  Fig.  5,  and  c  =  Ax/Ats. 
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FIG.  1.  Schematic  diagram  of  active  wall. 

FIG.  2.  Longitudinal  extension  of  piezoelectric  element. 

FIG.  3.  Typical  active-wall  cell  configuration  (not  to  scale). 

FIG.  4.  Isometric  drawing  of  active-wall  assembly  (not  to  scale). 

FIG.  5.  Block  diagram  of  master  control  circuit. 

FIG.  6.  Block  diagram  of  typical  channel  circuit. 

FIG.  7.  Comparison  of  load  line  and  safe  operating  areas  for  bipolar 
and  VMOSFET  power  transistors. 

FIG.  8.  Simplified  schematic  diagram  of  bridge  power  amplifier. 

FIG.  9.  Wave  amplitude  calibration  of  active  wall  for  a  frequency  of 
40  Hz. 

FIG.  10.  Wave  amplitude  response  of  active  wall:  O  100  V  input,  10.5um 
amplitude  at  10  Hz;  □  50  V  input,  4.2  urn  amplitude  at  10  Hz. 

FIG.  11.  Phase  response  of  active  wall  at  f  =  40  Hz,  Ats  *  0.76  ms, 
ax  *  1.25  mm,  c  *  166  cm/s,  and  X  *  4.1  cm:  straight  line,  control 

settings;  O  electro-optical  displacement  measurements. 


lBLE  IV.  Design  criteria  for  active  wall 


Quantity 

Symbol 

Value 

Units 

aximum  frequency 

f 

150 

Hz 

Lnimum  wavelength 

A 

10 

mm 

aximum  wave  speed 

c 

1.5 

m/s 

aximum  amplitude 

a 

25 

pirn 

ctive  length 

L 

40 

mm 

ell  spacing 

Ax 

1.25 

mm 

umber  of  cells 

L/Ax 

32 

umber  of  cells/wavelength  A /Ax 
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TABLE  III.  Physical  dimensions  of  membrane  and  Tolmien-Schlichting  waves  in 
a  water  laminar  boundary  layer. 


Source 

6 

(mm) 

X 

(cm) 

0  CD 

(cm/s) 

f 

(Hz) 

X 

(cm) 

c 

(cm/s) 

Non-Parallel3 

2 

7.7 

70 

27 

1.1 

30 

Parallel3 

2 

10.1 

91 

30 

1.2 

36 

Membrane 

2 

27.6 

251 

100 

.  1.3 

126 

Experiment3’15 

4.1 

14.5 

32 

6.5 

2.0  . 

13 

b  c 

Experiment 

1.5 

14.0 

234 

95 

0.9 

83 

a  At  Critical  Reynolds  Number 
b  3  r 

Keference  20 

c  At  Re5  =  3720  and  maximum  instability 


TABLE  II.  Nondimensional  characteristics  of  membrane  and  Tolmien-Schl ichting 
waves  in  a  laminar  boundary  layer. 


Source 

6U./V 

a  £ 

c/U® 

Non-Parallel  Theory3 

1400 

0.48 

1.15 

0.42 

Parallel  Theory^ 

1820 

0.42 

1.05 

0.40 

r 

Membrane  Model 

5000 

0.5 

1.0 

0.5 

a  Reference  16  at  critical  Reynolds  number 
b  Reference  14  at  critical  Reynolds  number 
c  Reference  17  at  reference  Reynolds  number 


a  X  L  ax  c  f  vmax 

Author  (mm)  (mm)  (cm)  (mm)  (m/s)  (Hz)  (cm/s) 


Kendall3 

3.18 

102 

122 

25.4 

3 

30 

59 

Weinstein  & 
Balasubramanianb 

0.02 

1.8 

40 

•  •  • 

0 

2,000 

25 

Wehrmannc 

0.001 

42 

4.2 

4.2 

2.1 

50 

0 

Present  Work 

0.013 

10 

4 

1.25 

1.2 

120 

1 
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Recent  numerical  calculations  by  Kuhn  et.  al.,27  indicate 
that  drag  reduction  may  be  possible  in  a  turbulent  boundary  with 
small  traveling  surface  waves.  However,  these  calculations  imply 
that  wavelength  must  be  less  than  1  mm  in  water  with  X/a  =  40. 

With  the  materials  in  the  present  active  wall,  wavelength  can 
be  reduced  to  approximately  3  mm  by  removal  of  one  insulator  sheet 
from  each  cell  and  reduction  in  the  number  elements  per  wavelength 
to  four  where  the  number  of  elements  per  wavelength  is  controlled 
electronically  through  selection  of  wave  speed  and  frequency. 

With  application  of  multilayered  monolithic  ceramic  technology, 
wavelength  could  be  reduced  to  0.25  mm,  but  such  technology  is 
not  yet  commercially  available.  When  multilayer  devices  become 
available,  a  single  piezoelectric  element  could  consist  of  ten 
layers  of  ceramic  separated  by  metal  electrodes  on  the  order  of 
a  few  molecules  thick.  Since  the  thickness  in  Eq.  (10)  would  be 
reduced  by  a  factor  of  ten,  voltage  could  be  reduced  by  an  equal 
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factor  for  the  same  amplitude. 

Significant  advantages  of  the  present  device  within  its  limi¬ 


tations  are  that  amplitude,  wave  speed,  and  frequency  can  be  con- 
tolled  independently.  Additionally,  the  wave  forms  are  not  limited 


to  sinusoids.  Since  the  drive  system  is  electronic,  a  traveling 
surface  wave  of  essentially  arbitrary  shape  can  be  produced  so 
that  compliant-wall  or  active-wall  theory  based  on  nonsinusoidal 


surface  waves  could  be  tested  with  the  device  described  in  this 
paper . 
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Phase  shift  measured  on  the  active  wall  by  the  displacement  sensor  is 


<J>  =  360  if  (12) 

where  t  is  the  measured  phase  delay  in  seconds  at  station  x  rela¬ 
tive  to  x  =  0 .  As  Fig.  11  indicates,  measured  phase  shift  from 
Eq.  (12)  is  in  agreement  with  design  values  determined  from  Eq.  (11). 
IV.  CONCLUSIONS 

An  active  wall  which  generates  short-wavelength  traveling 
surface  waves  has  been  successfully  designed  and  tested.  The 
active  wall  essentially  meets  the  design  requirements  outlined 
in  Table  IV,  although  the  maximum  displacement  amplitude  was  approx¬ 
imately  half  of  the  design  value. 

The  electronic  drive  produced  a  300-V  amplitude  signal  across 
a  4.0-kft  resistive  load  or  0.44-pF  capacitor.  For  the  resistive 
load  the  power  amplifier  had  flat  frequency  response;  however, 
with  a  capacitive  load  response  was  identical  with  that  of  a  current- 
fed  parallel  resonant  circuit  with  a  center  frequency  of  30.5  Hz 
and  a  quality  factor  (Q)  of  2.56.  In  the  present  work,  this  anomaly 
was  not  a  problem,  but  it  should  be  corrected  in  future  applications. 
Although  response  of  the  amplifier  circuit  driving  piezoelectric 
ceramics  was  similar  to  that  observed  with  a  capacitive  load, 
distortion  caused  by  saturation  of  the  ceramic  elements  appeared 
in  the  signal  to  the  elements  between  80  and  90- V  amplitude.  A 
possible  solution  to  this  problem  may  be  in  selection  of  a  piezo¬ 
electric  ceramic  with  a  lower  d3]_  constant  but  with  more  favorable 
values  of  large-signal  characteristics. 
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Design  of  a  Low-Turbulence  Water  Tunnel 

By 

Joel  T.  Park 

Southwest  Research  Institute 
San  Antonio,  Texas  78284 

INTRODUCTION 

A  low-turbulence  water  tunnel  was  designed  and  constructed  for  high 
lit y  fluid  mechanics  research.  A  test  section  was  constructed  for  boundary 
er  measurements  with  an  active  wall  for  both  laminar  and  turbulent  flows, 
facility  design  incorporated  an  available  pump  and  hydraulic  drive  motor, 
lure  1  is  a  scale  drawing  of  the  water  tunnel  while  Figure  2  is  a  photograph 
the  final  assembly.  The  nomenclature  of  the  various  components  is  identified 
Figure  1.  The  test  section  was  nominally  designed  for  a  maximum  velocity  of 
m/s  (32.8  ft/s),  and  its  cross  section  is  rectangular  with  dimensions  of 
2  x  22.9  x  76.2  cm  (6  x  9  x  30  in). 

The  design  features  of  a  number  of  facilities  were  reviewed  for  appli- 
lility  to  the  new  water  tunnel.  Some  typical  tunnels  are  listed  in  Table  1. 
ly  of  the  features  of  the  Case  tunnel,  where  a  number  of  successful  laminar 
indary  layer  experiments  have  been  performed  [3,4],  and  the  smaller  Penn  State 
inel  [7]  were  included  in  the  design.  From  Table  1,  a  reasonable  design  goal 
>eared  to  be  a  relative  turbulence  intensity  of  0.1%  in  the  longitudinal  com- 
lent.  Additionally,  some  of  the  more  recent  concepts  in  contraction  and  set- 
ng  chamber  design  have  been  incorporated  to  ensure  optimum  performance. 

Since  the  SwRI  tunnel  is  relatively  small,  the  cost  of  the  materials 
not  a  major  factor  in  the  design.  The  Case  tunnel  [2],  which  is  about  the 
le  size  as  the  SwRI  tunnel,  is  constructed  from  all  aluminum.  For  the  same 
'ength  in  the  design,  thinner  304  stainless  plate  is  required;  consequently, 

>  total  cost  of  304  stainless  steel  and  aluminum  is  nearly  the  same.  Also, 
linless  steel  is  more  corrosion  resistant;  therefore,  it  was  selected  as  the 
mary  material  for  the  tunnel.  The  design  pressure  for  the  facility  is  207  kPa 
I  psi). 
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FIGURE  1.  SCALE  DRAWING  OF  LOW- TURBULENCE  WATER  TUNNEL  AT  SOUTHWEST  RESEARCH  INSTITUTE 


' 


PHOTOGRAPH  OF  LOW-TURBULEN 


TABLE  5.  INFLUENCE  OF  CONTRACTION  ON  FREESTREAM  TURBULENCE 
FOR  A  CONTRACTION  AREA  RATIO  OF  16 


Hussain  [25]  reported  that  exit  relative  turbulence  intensity,  (ou/U)e,  was 
pendent  of  the  inlet  conditions  from  grid  turbulence.  Conclusions  from  the 
recent  results  of  Tan-atichat  [28]  and  Nagib,  et  al.  [29]  which  were  pri- 
ly  from  matched-cubic  nozzles  include  the  following: 


Optimum  contraction  ratio  is  9 


L/D-j  should  be  one 


Upstream  integral  scale  and  (ou/U)i  influence 
(ou/u)e  and  other  exit  turbulence  characteristics 
for  large  contractions 


d.  Differences  between  matched-cubic  and  fifth  order  contours 
are  insignificant. 

lough  some  consensus  exists  among  various  researchers,  some  controversy  remains 
the  effects  of  contractions  on  turbulence. 


Nagib,  et  al.  [29]  have  formulated  design  charts  for  optimum  turbulence 
ipulation  with  a  matched-cubic  nozzle  of  contraction  ratio  9.  These  charts  are 
ited  to  longitudinal  turbulence  intensities  of  0.14  <  (ou/U)e  <  0.32  at  the  noz- 
exit.  Although  these  charts  overlap  the  range  of  interest,  0.1  <  (ou/U)e  <  0.2, 
>  information  became  available  after  the  contraction  was  designed  and  fabricated. 


'j 


On  the  basis  of  the  results  of  these  previous  investigators,  a  matched- 
ic  profile  was  selected  for  the  present  application.  The  contraction  ratio  was 
»en  as  16.67.  This  choice  allows  the  flow  within  the  settling  chamber  to  be 
te  the  critical  Reynolds  number  for  the  screens  and  to  maintain  the  loading  on 
screens  within  their  structural  limits.  The  turbulence  characteristics  for 
ious  nozzles  with  C  =  16  are  summarized  in  Table  5.  Tan-atichat  [28]  provides 
following  correlation  for  the  turbulence  reduction  factor  of  a  matched-cubic 
File 

°ue/°ui  -  C(C-5)/6]*  (12) 


6  <  C  <  30. 
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b.  6  screen  frames  including  an  entrance  screen  for  the  back 
pressure  on  the  wide-angle  diffuser 

c.  Blank  section  22.9  cm  long  (72  honeycomb  mesh  lengths) 

The  blank  section  and  two  screen  frames,  one  adjacent  to  the  honeycomb 
le  second  downstream  of  the  blank  section,  provide  a  decay  distance  of 
=  120  for  the  turbulence  between  the  honeycomb- screen  combination  and  the 
downstream  screen.  Also,  each  of  the  internal  components  is  separated 
3.2  mm  (0.125  in)  rubber  gasket. 

©  Contraction  Section 

A  contraction  section  was  designed  in  conjunction  with  the  settling 
er  for  minimization  of  the  turbulence  in  the  test  section.  Linear  rapid 
rtion  theory  [23]  indicates  that  the  turbulence  reduction  factor  decreases 
onically  with  an  increase  in  contraction  ratio.  For  large  contraction  ratios 
esult  is 


(oue/oui)'  *  [3/(4CJ)]tln(4C3)-l]  (11) 

i  C  is  the  contraction  area  ratio.  Since  turbulence  reduction  factor  goes 
iro  for  large  contraction  ratios,  eq.  (11)  implies  that  large  contraction 
is  are  required  for  low-turbulence  facilities. 

Experimentally,  Uberoi  [24]  and  others  demonstrated  that  the  theory  in 
11)  is  incorrect.  The  experimental  results  of  Ramjee  and  Hussain  [25] 

:ate  that  the  optimum  contraction  ratio  is  40  with  no  significant  reductions 
irbulence  above  about  20.  Tsuge  [26]  has  recently  published  a  theory  which 
i  agreement  with  experiment.  His  theory  accounts  for  the  amplification  of 
t  eddies  which  are  ignored  by  the  linear  theory. 

In  addition  to  contraction  ratio,  other  key  factors  in  nozzle  design  are 
i,  length  to  inlet  diameter  ratio  ( L/D ^ ) ,  and  the  characteristics  of  the  up- 
im  turbulence.  In  their  experiments,  Hussain  and  Ramjee  [27]  concluded  that 
latched-cubic  contour  was  the  best  of  four  nozzles  they  tested,  and  Ramjee 
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>teel  screen  that  was  purchased  has  a  mesh  size  of  0.67  cm  (1/38  in)  and  a  wire 
Jiameter  of  0.114  mm  (4.5  mil).  Thus,  6  =  68.7%,  Re<j  =  68.3,  and  Re<jM  *  99.4, 
»ihile  the  critical  Rej  from  Schubauer,  et  al .  [21]  is  57.2  (Re^  =  83.2). 

The  reduction  in  turbulence  for  the  flow  through  n  screens  is  given  by 
the  turbulence  reduction  factor  from  Dryden  and  Schubauer  [22]  as 


aun/aui  =  (1  +  k)“n/2 


(9) 


while  Loehrke  and  Nagib  [17]  have  observed  the  correlation 

<W°ui  =  (1  +  k)‘n/2-7.  (10) 

Equation  (4)  provides  an  estimate  of  the  pressure  drop  coefficient.  The  estimates 
of  k  from  this  correlation  and  others  are  listed  in  Table  4.  The  turbulence  re¬ 
duction  factors  from  eqs.  (9)  and  (10)  for  k  =  0.662  and  four  screens  (n  =  4)  are, 
respectively,  0.362  and  0.471.  On  the  basis  of  these  estimates,  the  total  tur¬ 
bulence  reduction  factor  for  the  settling  chamber  is  0.061  to  0.15. 

The  screen  frames  were  fabricated  from  304  stainless  steel  square  tube 
with  sides  of  7.62  cm  (3  in)  and  wall  thickness  of  3.2  mm  (0.125  in).  The  frames 
have  sufficient  strength  to  support  the  loading  of  the  screens,  and  they  provide 
a  space  of  l/W  =  114  for  the  decay  of  turbulence  between  screens.  The  loading 
of  the  screens  was  estimated  from  Domholdt  [2].  The  screens  are  preloaded  by 
1750  N/m  (10  lb/in)  in  one  direction  and  bonded  to  the  frames.  With  this  pre¬ 
load,  the  screen  deflection  is  estimated  to  be  less  than  one  percent  of  the 
channel  width. 

The  settling  chamber  is  a  square  shell  of  6.4  mm  (0.25  in)  stainless  plate 
with  flanges  and  ribs.  The  internal  components  are  sandwiched  between  the  flanges 
of  the  wide-angle  diffuser  and  contraction  section.  The  flow  through  area  is 
76.2  cm  square  by  76.2  cm  (30  in)  long.  The  internal  components  include  the 
following: 


a.  Honeycomb  with  hexagonal  cells,  3.2  mm  x  5.08  mm  long 


installed.  The  honeycomb  eventually  corroded  away,  and  it  was  replaced  by 
of  stainless  steel  [6].  The  honeycomb  at  Penn  State  has  hexagonal  cells  with 
idth  across  flats  of  0.56  cm  by  47.62  cm  long  (0.22  x  18.75  in)  or  an  £/M  of  85. 
cells  in  the  NBS  tunnel  [5],  which  are  also  stainless  steel,  are  0.32  x  61  cm 
125  x  24  in)  or  an  1/ M  of  192.  The  theoretical  requirement  for  the  large  i/M 
i]  is  based  upon  fully  developed  turbulent  flow  in  the  cells. 

Since  the  cost  of  thick  honeycomb  appeared  to  be  excessive  for  the  present 
>lication,  a  more  economical  approach  with  a  combination  of  screens  and  honey- 
ob  suggested  by  Loehrke  and  Nagib  [17,19]  was  adopted.  This  method  has  been 
plied  in  modifications  to  the  Case  tunnel  which  now  contains  screens,  honey- 
mb,  and  foam  [4]. 

The  important  design  parameters  in  the  selection  of  honeycomb  are  length 
mesh  size  ratio  (£/M)  and  the  mesh  Reynolds  number,  ReM  =  UsM/v,  where  Us  is 
e  settling  chamber  velocity.  Honeycomb  is  effective  in  removing  swirl  and 
mpening  the  lateral  turbulence  scales,  and  it  should  have  l/M  >  10  [19].  The 
at  exchanger  functions  as  a  honeycomb  for  the  removal  of  large  scales  while 
e  honeycomb  in  the  settling  chamber  manipulates  primarily  the  smaller  scales, 
e  honeycomb  in  the  settling  chamber  is  fabricated  from  304  stainless  steel  foil 
th  a  thickness  of  76  pm  (3  mil),  and  the  cell  size  is  3.2  mm  by  50.8  mm  long 
.125  x  2  in)  1/M  =  16.  The  Reynolds  number,  Re#,  at  the  design  test  section 
locity  is  190o. 

One  of  the  honeycomb  configurations  (plastic  soda  straws)  tested  by 
ehrke  and  Nagib  [17,19]  was  near  the  design  condition  of  the  present  honey- 
mb.  In  those  experiments  a  screen  near  the  honeycomb  at  Ax/M  <5.7  signifi- 
ntly  reduced  the  standard  deviation  or  rms  velocity,  oue/ou-j,  or  turbulence 
duction  fa  tor  was  0.31  without  the  screen  and  0.17  with  the  screen. 

According  to  Lumley  and  McMahon  [18],  the  decay  rate  of  turbulence  through 
reens  and  honeycomb  at  a  turbulent  Reynolds  number  is  higher.  A  recent  review 
flow  through  screens  is  by  Laws  and  Livesey  [20].  The  critical  Reynolds  number 
r  screens  is  based  upon  wire  diameter,  d.  For  a  Reynolds  number  based  upon  grid 
t  velocity,  the  critical  value  of  Re<jM  =  Usd/(Bv)  is  80  [12]  while  Schubauer, 
al.  [21]  measured  critical  Re^  =  Usd/\;  as  a  function  of  solidity.  The  stainless 
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23.9°C  (75°F).  The  estimated  power  dissipation  of  the  heat  exchanger  is  10.6  kW 
(36,000  Btu/hr). 

The  heat  exchanger  is  constructed  from  all  aluminum.  The  tube  bundle 
consists  of  247  tubes  with  an  outside  diameter  of  2.54  cm  (1  in),  wall  thickness 
of  0.89  mm  (0.035  in),  and  length  of  40.6  cm  (16  in).  The  tubes  are  welded  to 
19-mm  (0.75  in)  thick  tube  sheets  which  also  serve  as  flanges.  The  hole  pat¬ 
tern  in  the  tube  sheet  is  a  equilateral  triangle  with  center  to  center  distance 
of  3.18  cm  (1.25  in).  The  shell  is  also  19-mm  (0.75  in)  plate  with  an  inside 
square  cross  section  of  48.9  cm  (19.25  in). 

As  an  effective  turbulence  manipulator,  the  heat  exchanger  should  have 
the  largest  possible  open-area  ratio,  6.  According  to  Morgan  [16],  the  open- 
area  ratio  of  screen  should  be  larger  than  0.5;  otherwise,  a  flow  instability 
will  occur.  Also,  Loehrke  and  Nagib  [17]  have  noted  anomalous  behavior  for 
low-  6(8  =  0.30)  perforated  plates.  From  geometry,  the  local  open-area  ratio 
for  a  perforated  plate  with  a  triangular  hole  pattern  is 

Bf  0rVl>6)(d/£)2  (7) 

where  d  is  the  hole  diameter  and  £  is  the  centerline  distance  between  holes, 
and  the  average  value  is 

6  =  mrd2/(4A)  (8) 

where  n  is  the  number  of  holes  and  A  is  the  inside  cross-sectional  area  of  the 
mating  ducts,  the  wide-angle  diffuser  and  large  turning  section.  The  local  and 
average  values  of  the  present  design  are,  respectively,  50%  and  42%. 

(D  Settling  Chamber 

More  recently  turbulence  in  water  tunnels  has  been  managed  with  instal¬ 
lation  of  honeycomb  in  the  settling  chamber.  The  large  water  tunnel  at  Penn 
State  required  large  wire-diameter  screens  on  the  basis  of  strength  requirements 
which  caused  vortex  shedding  problems  [18],  Consequently,  a  design  criteria  was 
developed  by  Lumley  and  McMahon  [18]  for  honeycomb,  and  an  aluminum  honeycomb 


of  zero.  The  gap/chord  ratio  is  slightly  larger  for  thicker  vanes  [15].  Pope  [9] 
recommends  a  tunnel  width/gap,  w/h,  ratio  greater  than  six  and  a  large  chord.  For 
convenience,  the  large  turning  section  was  divided  into  twenty  increments  with  19 
vanes,  and  the  small  turning  section  into  twelve  with  11  vanes.  The  resulting 
increment  across  the  width  of  the  turning  sections  is  2.54  cm  (1  in).  The  chord 
was  then  selected  as  four  times  this  value  or  10.2  cm  (4  in).  From  geometry  the 
chord  and  radius  are  related  by 

c/r  =  2  sin  (  9  /2  )  (6) 

where  0  is  the  angle  of  the  arc  which  is  85°  for  an  optimum  turning  vane.  From 
eq.  (6)  the  radius  is  7.5  cm  (2.96  in). 

The  vanes  were  then  fabricated  from  15.2  cm  (6  in)  diameter  stainless-steel 
tubing  with  a  wall  thickness  of  3.2  mm  (0.125  in).  The  upper  and  lower  walls  of 
the  turning  sections  were  rolled  from  stainless-steel  plate  with  the  same  corner 
radius.  In  the  final  assembly,  the  gap/chord  ratio  across  the  corners  is  0.36. 

The  tunnel  width/gap  ratios  for  the  large  and  small  turns  are,  respectively,  14.1 
and  8.5.  The  plate  thickness  in  the  large  turn  Is  9.5  mm  (0.375  in),  and  in  the 
small  turn  6.4  mm  (0.25  in). 

(4)  Heat  Exchanger 

A  heat  exchanger  was  designed  for  removal  of  heat  from  frictional  heating. 
The  device  also  functions  as  a  turbulence  manipulator,  i.e.  coarse  honeycomb,  for 
dampening  large  scale  turbulence.  The  device  was  designed  with  approximately  the 
same  area  for  heat  transfer  as  that  of  Domholdt  [2]  with  some  additional  improve¬ 
ments  for  efficiency. 

The  heat  exchanger  essentially  functions  as  a  counter-flow  heat  exchanger. 
The  coolant  water  enters  a  manifold  at  the  bottom  downstream  end  of  the  heat 
exchanger,  is  routed  over  the  tube  bundle  by  two  vertical  baffle  plates,  and 
exits  a  second  manifold  at  the  top  upstream  end.  Preliminary  heat  transfer 
calculations  indicate  that  the  tunnel  temperature  will  be  maintained  at  26.6°C 
(80°F)  for  the  maximum  design  speed  with  a  coolant  flow  of  1.3  l/s  (20  gpm)  at 


50.8  cm  (20  in)  and  76.2  cm  (30  in),  or  the  area  ratio  is  2.25.  From  eq.  (3) 
the  number  of  screens  is  then  three  for  a  k  =  0.725.  The  dimensions  of  the 
first  two  screens  were  computed  from  eq.  (2)  as  58.21  cm  (22.92  in)  and 
66.71  cm  (26.27  in). 

The  total  length  of  the  wide-angle  diffuser  was  chosen  as  71.8  cm  (28.25  in). 
The  first  screen  is  located  at  37.5  cm  (14.8  in)  from  the  entrance.  This  distance, 
which  is  approximately  16  tube  diameters  from  the  heat  exchanger,  was  selected  to 
ensure  adequate  mixing  of  the  jets  from  the  heat  exchanger  before  the  first  screen 
is  encountered. 

The  diffuser  was  designed  so  that  the  first  two  screens  could  be  removed 
without  emptying  the  tunnel.  The  third  screen  is  located  in  the  entrance  of  the 
settling  chamber.  The  wall  thickness  of  the  diffuser  is  9.5  mm  (0.375  in). 

( 9)  Split  Diffuser.  The  diffuser  downstream  of  the  test  section  is  similar 
in  design  to  that  of  Domholdt  [2].  The  inlet  and  outlet  inside  dimensions  are, 
respectively,  15.2  x  22.9  cm  (6  x  9  in)  and  30.5  cm  (12  in)  square.  The  total 
angle  is  16°  between  the  upper  and  lower  surfaces  with  a  3.2  mm  (0.125  in) 
splitter  plate  in  the  center  while  total  angle  for  an  equivalent  conical  diffuser 
is  9.9°.  The  length  of  the  divergent  section  and  splitter  plate  is  54.2  cm 
(21.35  in).  A  short  parallel  section  is  included  upstream  of  the  splitter  plate 
so  that  the  total  length  of  the  diffuser  is  63.5  cm  (25  in).  The  wall  thickness 
is  6.4  mm  (0.25  in). 

dll)  Transition.  The  transition  in  the  vertical  return  has  the  same  inlet  and 
outlet  areas.  The  dimensions  of  the  inlet  and  outlet  are,  respectively,  30.5  cm 
(12  in)  square  and  34.3  cm  (13.5  in)  in  diameter.  The  section  is  108  cm  (42.5  in) 
long  with  a  6.4  mm  (0.25  in)  wall  thickness. 

(7)  4  (J§)  Turning  Sections 

According  to  Bradshaw  and  Pankhurst  [12]  and  Pankhurst  and  Holder  [15], 
the  most  efficient  turning  section  consists  of  turning  vanes  fabricated  from 
sheet  metal  with  a  circular  arc.  The  optimum  vane  has  a  gap/chord,  h/c,  ratio 
of  0.25  for  thin  vanes,  a  leading  edge  angle  of  5°,  and  a  trailing  edge  angle 
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The  wall  thickness  of  the  large  diffuser  is  9.5  mm  (0.375  in). 


0  Wide  Angle  Diffuser.  A  30°  wide-angle  diffuser  was  incorporated  in  the 
design  upstream  of  the  settling  chamber  for  a  reduction  in  length.  Screens  were 
included  to  prevent  flow  separation  and  to  reduce  the  level  of  turbulence.  The 
section  was  designed  by  the  method  of  Schubauer  and  Spangenberg  [10]  so  that  the 
pressure  recovery  was  balanced  by  the  pressure  loss  of  the  screen.  The  area  ratio 
at  the  screen  location  is  given  by  [10] 

Ai/Ai-i  *  (ki  +  1)*  (2) 

i  =  1,2, . n 

where  Aj  is  the  area  of  the  ith  screen  and  ki  is  the  pressure  drop  coefficient 
of  the  screen.  For  n  screens  with  the  same  ki  the  total  area  ratio  of  the 
diffuser  is 


An/A0  •  (k  +  l)n/2.  (3) 

The  pressure  drop  coefficient  recommended  by  Nagib  [11]  is  given  by 

k  =  (1  -  S )/ B2  (4) 

which  is  an  unpublished  result  valid  for  0.5  <8<  0.7  where  the  open-area  ratio 
of  the  screen  is 


6  •  (1  -  d/M)2  (5) 

d  is  the  wire  diameter  and  M  is  the  mesh  length.  The  open-area  ratio  is  related 
to  solidity,  o,  by  B*  1  -  o.  Other  available  correlations  for  the  pressure  drop 
correlation  are  summarized  in  Table  4,  but  screens  will  be  discussed  in  more  detail 
in  a  subsequent  section  on  the  settling  chamber. 

For  the  purpose  of  the  design  calculations,  a  mesh  size  of  0.635  mm  (1/40  in) 
and  a  wire  diameter  of  0.114  mm  (0.0045  in)  were  selected  although  a  slightly 
larger  mesh  was  purchased.  The  inlet  and  outlet  dimensions  of  the  diffuser  are 
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TABLE  2.  PERFORMANCE  CHARACTERISTICS  OF 

BERKELEY  VERTICAL  MIXED  FLOW  PUMP 
MODEL  M12-I2 


Flowrate 

0.315 

m3/s 

(5000  gpm) 

Rotational  Rate 

1760 

rpm 

Power  Required 

50 

kW 

(67  hp) 

Submergence  Required 

0.6 

m 

(2  ft) 

Total  Dynamic  Head 

125 

kPa 

(42  ft  of  water) 

TABLE  3.  PERFORMANCE  CHARACTERISTICS 
DENISON  VANE  TYPE  HYDRAULIC 
MODEL  M3D-138-21N 

OF 

MOTOR 

Maximum  Power 

65 

kW 

(87  hp) 

Rotational  Rate 

2000 

rpm 

FI owrate 

*  4.7 

*/s 

(75  gpm) 

Maximum  Pressure 

13.8 

MPa 

(2000  psi) 
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DESIGN  DETAILS 


©  Pump1 

The  water  tunnel  was  designed  around  an  available  Berkeley  mixed  flow 
pump.  The  pump  is  powered  with  a  Denison  hydraulic  motor  so  that  tunnel  speed 
is  controlled  by  hydraulics.  Speed  control  can  be  accomplished  with  a  manual 
valve  or  with  a  Koehring-Pegasus  model  1282A  two-stage  servovalve  for  remote 
control.  The  specifications  of  the  pump  and  motor  are  summarized  in  Tables  2 
and  3.  The  total  dynamic  head  of  the  pump  is  more  than  5  times  the  estimated 
pressure  loss  of  the  tunnel. 

The  pump  is  fitted  with  American  Standards  Association  (ASA  B16.1  and 
B16.5)  slip-on  flanges  for  14-inch  diameter  pipe  (35.6  cm).  The  return  leg 
of  the  water  tunnel  at  pump  inlet  contains  a  14-inch  diameter  pipe  ©  and 
standard  long-radius  elbow  ©  .  The  elbow,  pipe,  and  pump  are  galvanized  steel. 

Di f fusers 

©Large  Diffuser.  The  optimum  divergence  angle  for  minimum  power  loss  in 
diffuser  is  given  by  Pope  [9]  as 

tan  (a/2)  =Vf/4.8  (1) 

where  a  is  the  total  divergence  angle  and  f  is  the  pipe  friction  factor.  In  the 
present  application,  the  optimum  divergence  angle  is  about  6°.  In  the  opinion 
of  Dumholdt  [2],  an  8°  angle  is  adequate;  consequently,  the  8°  value  was  accepted 
to  reduce  the  height  of  the  tunnel. 

The  diffuser  also  serves  as  a  transition  section  from  a  circle  to  a 
square.  The  inlet  and  outlet  inside  dimensions  are,  respectively,  34.2  cm 
(13.5  in)  diameter  and  50.8  cm  (20  in)  square.  The  resulting  length  for  an  8° 
divergence  of  the  sides  is  1.18  m  (46.5  in). 

Circled  numbers  correspond  to  those  of  the  tunnel  drawing  in  Figure  1. 


As  Figure  2  indicates,  the  water  tunnel  is  supported  vertically  by  three 
large  steel  frames,  a  small  steel  frame,  and  the  support  stand  for  the  pump. 

The  dry  mass  of  the  facility  excluding  the  support  structure  is  2,400  kg 
(5,300  lbs),  the  total  mass  excluding  the  support  structure  is  4,900  kg  (10,700 
lbs),  and  the  water  volume  is  2.45  m3 (650  gal)  .  The  tunnel  is  supported  on 
twenty  Uni  sorb  model  LR750  machinery  mounts  with  leveling  screws  for  vibration 
isolation.  The  total  height  of  the  tunnel  centerline  to  centerline  is  1.9  m 
(73.5  in)  while  the  horizontal  centerline  length  is  5.2  m  (204  in). 


20 

25 

a> 

CT»  CT> 

7.5 

&*-  *  ' 

'  -• 

OV 

oo 

£■*  vo 

iH 

V 

ID 

o» 

X 

o 

*3* 

CO 

X 

X 

CO 

X 

O 

CVJ 

X 

_  Q 

o 

CO 

!_  0  ^ 

VO 

X 

o 

ID 

VO 

• 

CVJ  . 

.  - ".  *  .  ' 

ir> 

ir> 

fH 

CVJ  o 

o 

CO 

•-H 

VO 

r-l  ro 

CO 

From  the  range  of  experimental  values  in  Table  5,  the  estimated  turbu¬ 
lence  reduction  factor  for  the  contraction  and  settling  chamber  combination  is 
0.065  <  Oue/°ui  £  0-34  or  the  ratio  of  relative  intensities  is  0.0039  <  (c\|/U)e/ 
(ou/U)i  <  0.020.  Thus  if  (ou/U)j  entering  the  settling  chamber  is  no  more  than 
10%,  the  relative  turbulence  intensity  of  the  water  tunnel  should  be  within  the 
design  goal  of  0.2%  in  the  test  section  for  a  properly  designed  contraction. 

The  contraction  was  designed  from  the  procedure  of  Morel  [31,32].  Morel's 
design  curves  for  an  axi symmetric  nozzle  [31]  were  from  a  numerical  finite- 
difference  solution  of  Euler's  equations.  He  performed  calculations  on  six 
power-law  shapes,  orders  two  through  five,  and  concluded  that  the  matched-cubic 
shape  was  the  best.  The  design  procedures  have  the  following  restrictions: 

Cpi  >  0.1  0.75  <  L/Di  <  1.25 

Cpe  <0.1  0.2  Xm/L  ^  0.8 

2  <  C  <  25 

where  L  is  the  contraction  length,  Di  is  the  inlet  diameter,  and  Xm  is  the 
match  point  of  cubic  curves.  The  pressure  coefficients  Cpe  and  Cpi-  are  defined 
as  follows: 


pe  =  1  -  (Uoo/Uy^g) 

(13a) 

pi  =  1  "  (Uwi /^s) 

(13b) 

where  Uoois  the  test  section  velocity,  Us  is  the  settling  chamber  velocity,  and 
Uwi  and  Uwe  are  the  inlet  and  exit  velocities  at  the  wall.  The  pressure  coeffi¬ 
cients  are  a  measure  of  the  velocity  profile  flatness,  and  they  are  chosen  in 
the  design  method. 


The  following  is  a  brief  summary  of  the  design  method  where  equation  and 
figure  numbers  refer  to  Morel  [31]: 


a.  Choose  C,  Cpe,  and  Cp-j 


b.  Compute  equivalent  diameters  from  D  =  2  Va/tt 

c.  Diminish  Cp-j  by  302 

d.  Determine  Fe  and  Gj  from  Fig.  10  and  11 

e.  Compute  Xm/L  from  eq.  (10)  or  Fig.  12 

f.  Compute  L/Di  from  eq.  (8b) 

g.  Check  pressure  gradient  for  possible  separation 

h.  Add  straight  sections  to  the  entrance  of  at  least 
0 . 2D-j  and  to  the  exit  of  0.3De. 

The  total  nozzle  length  is  then 


Lt  =  L  +  a  i  D-j  +  agDg.  (14) 

Since  the  velocity  is  relatively  low  in  this  application,  the  critical 
separation  criterion  may  be  the  influence  of  an  adverse  pressure  gradient  on  a 
laminar  boundary  layer.  The  boundary  layer  separation  criterion  is  from  Stratford 
[33],  and  pressure  coefficient  form  for  laminar  separation  in  contraction  coordi¬ 
nates  is  given  by  eq.  (14)  in  Morel  [32].  Separation  criteria  determined  the 
entrance  length  between  the  contraction  and  last  screen  in  the  present  design. 

The  equations  for  the  contraction  coordinates  are  given  by  eqs.  (10a)  and 
(10b)  in  Morel  [32].  The  versions  of  these  equations  in  Morel  [31],  eqs.  (5a)  and 
(5b),  are  incorrect.  The  dimensions  of  the  final  design  are  as  follows: 


Equivalent  inlet  diameter 

85.98  cm 

33.85 

Equivalent  exit  diameter 

21.06 

8.29 

Match  point 

74.52 

29.34 

Contraction  Length 

103.75 

40.85 

Entrance  Length 

17.78 

7.0 

81 


Exit  length 
Total  length 


6.35 

127.86 


2.5 

50.34 


The  contraction  section  was  rolled  from  6.4  mm  (0.25  in)  stainless 
steel  plate  and  welded.  The  profile  and  layout  coordinates  for  fabrication 
were  generated  by  a  computer  program.  Obtaining  the  layout  coordinates  for  the 
cubic  curve  required  a  numerical  integration.  The  interior  surface  of  the  con¬ 
traction  was  polished  to  a  measured  finish  of  0.30  pm  (12  p  in)  rms. 

(8)  Test  Section 

The  cross-sectional  area  of  the  tunnel  was  determined  by  the  design  speed 
of  10  m/s  and  maximum  flowrate  of  the  pump.  The  dimensions  are  15.2  cm  high  by 
22.9  cm  wide  by  76.2  cm  long  (6  x  9  x  30  in).  The  aspect  ratio  of  the  cross- 
section  is  near  the  VF value  for  traditional  wind  tunnels  [12],  and  the  length 
is  sufficient  for  most  boundary  layer  research.  The  tunnel  design  will  readily 
accommodate  other  test  section  lengths.  Figure  3  is  a  photograph  of  the  assembled 
test  section. 

The  test  section  includes  a  flat  plate  with  a  bleed  slot  for  removal  of 
the  tunnel  boundary  layer.  The  concept  is  similar  to  that  of  the  Case  tunnel  where 
several  successful  boundary  layer  experiments  have  been  performed  [3,4].  This 
design  was  chosen  because  it  will  allow  installation  of  externally  bulky  devices 
such  as  active-wall  devices,  manifolds  for  boundary-layer  injection,  surface- 
mounted  hot-film  probes,  and  electro-magnets  for  driving  vibrating  ribbons. 

The  top  and  sides  of  the  test  section  are  fabricated  from  Plexiglas  for 
flow  visualization  and  laser  Doppler  anemometer  (LDA)  measurements.  The  bottom 
flat  plate  was  machined  from  aluminum  and  plated  with  Tufram  L-4,  a  teflon  coat¬ 
ing  from  General  Magnaplate  Corporation.  The  flow  side  of  the  plate  was  polished 
to  a  finish  of  0.15  pm  (6  p  in)  rms,  but  the  final  finish  with  the  Tufram  is 
0.38  pm  (15  p  in)  rms.  The  nose  of  the  plate  is  an  ellipse  with  a  length/ 
thickness  ratio  of  2.86,  and  it  is  located  5.08  cm  (2  in)  from  the  entrance. 

The  dimensions  of  the  plate  are  69.2  cm  long  by  22.9  cm  wide  by  1.78  cm  thick 
(27.25  x  9  x  0.7  in).  The  flanges  were  machined  from  1.9  cm  (0.75  in)  aluminum 


PHOTOGRAPH  OF  TEST  SECTION  IN  LOW- TURBULENCE 
WATER  TUNNEL  AT  SOUTHWEST  RESEARCH  INSTITUTE 


plate  and  anodized.  Since  the  flat  plate  protrudes  into  the  flow,  the  exit 
dimensions  of  the  test  section  are  13.7  cm  high  and  22.9  cm  wide  (5.4  x  9  in) 
All  mating  surfaces  are  sealed  with  0-rings  so  that  the  test  section  can  be 
readily  disassembled. 


CONCLUSIONS  AND  RECOMMENDATIONS 


A  low-turbulence  water  tunnel  has  been  designed  and  constructed  for 
boundary  layer  experiments.  The  latest  available  information  on  design  has 
been  incorporated  in  the  facility.  The  nominal  design  velocity  is  10  m/s 
(32.8  ft/s).  The  test  section  has  three  Plexiglas  sides  for  flow  visualization 
and  LDA  measurements,  and  the  fourth  side  is  a  polished  aluminum  plate  with  a 
Tufram  coating  for  boundary  layer  measurements.  A  bleed  slot  is  included  in 
the  test  section  for  removal  of  the  tunnel  boundary  layer. 

Turbulence  in  the  test  section  and  its  flow  profile  are  controlled  by 
a  settling  chamber  and  matched-cubic  contraction  section.  The  settling  cham¬ 
ber  contains  six  stainless  steel  screens  with  a  fine  mesh  and  large  open  area 
ratio,  and  the  contraction  has  a  polished  interior  surface.  Estimates  on  the 
basis  of  other  experimental  results  indicate  that  the  relative  turbulence  inten¬ 
sity  in  the  test  section  should  be  between  0.04  and  0.2%  if  the  entrance  turbulence 
level  to  the  settling  chamber  is  less  than  10%.  According  to  Nagib  [11]  a  rela¬ 
tive  longitudinal  turbulence  intensity  of  0.2%  should  be  readily  attainable  with 
the  current  design.  Performance  characteristics  in  the  near  future  will  be 
measured  with  hot-film  anemometry. 

An  alternate  approach  in  tunnel  design  would  be  to  select  a  smaller 
contraction  ratio  and  thicker  honeycomb.  As  Table  1  indicates  other  successful 
facilities  have  been  built  with  a  contraction  ratio  of  9.  Although  the  cost  of 
the  honeycomb  would  increase  substantially,  potential  reductions  in  cost  would 
result  from  the  following: 

a.  Smaller  contraction 

b.  Smaller  settling  chamber 

c.  No  wide-angle  diffuser 

d.  No  screens. 

Preliminary  calculations  indicate  that  a  properly  designed  honeycomb 
without  screens  will  produce  the  same  relative  turbulence  intensity  in  the  test 


section  if  that  which  enters  the  settling  chamber  were  the  same.  Estimates 
for  the  current  design  with  C  !  9  indicate  that  a  cell  size  of  3.2  mm  by  152  mm 
( A/M  =  48)  would  be  reasonable.  The  Reynolds  numbers  based  on  mesh  size  and  length 
are  ReM  *  3,500  and  Reg,  =  1.7  x  10^.  From  Fig.  2  of  Lumley  and  McMahon  [18]  the 
pressure  drop  coefficient  is  2.2,  and  from  Fig.  1  the  turbulence  reduction  factor 
is  probably  less  than  0.2.  From  eq.  (12)  the  turbulence  reduction  factor  for  the 
contraction  with  C  a  9  is  0.82,  or  the  total  for  the  honeycomb  and  contraction 
will  be  less  than  0.16.  If  the  turbulence  level  entering  the  honeycomb  were  less 
than  10%,  the  relative  turbulence  intensity  in  the  test  section  would  be  less 
than  0.2%,  which  is  the  maximum  estimate  for  the  present  design  with  C  *  16.7. 

A  properly  designed  settling  chamber  requires  knowledge  of  the  upstream 
turbulence  characteristics  which  are  not  known  for  a  new  facility  design.  Also, 
the  results  from  various  experimenters  which  were  applied  in  the  current  calcu¬ 
lations  are  not  in  agreement.  For  example,  the  low  and  high  values  of  turbulence 
reduction  factor  for  the  various  components  yield  an  estimated  relative  turbulence 
intensity  between  0.04  and  0.2%  in  the  current  design  with  C  =  16.7.  Consequently, 
measurements  in  the  test  section  of  the  completed  facility  are  required  for  vali¬ 
dation  of  a  design. 

The  cost  of  the  honeycomb  for  a  C  =  9  design  with  a  cell  size  of  3.2  mm 
and  cell  length  of  152  mm  can  be  estimated  from  the  cost  of  the  honeycomb  for 
the  SwRI  tunnel.  At  65£  per  cubic  inch,  the  cost  for  the  C  *  9  design  would 
be  approximately  $2000  in  comparison  to  $1200  for  the  SwRI  tunnel.  Although 
the  price  of  the  honeycomb  would  double,  other  cost  benefits  may  be  possible 
from  the  simplicity  in  design  for  a  contraction  ratio  of  9. 
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ABSTRACT 


The  Kelvin-Helmholtz  solution  for  flow  over  a  traveling  wavy  surface 
which  extends  to  infinity  upstream  and  downstream  is  extended  to  the  case  of 
flow  over  a  semi-infinite  wavy  plate.  The  analytical  solution  is  obtained  by 
conformal  mapping  and  integral  transforms.  The  leading  edge  modifies  the 
Kelvin  solution  to  form  a  secondary  flow  pattern  which  appears  to  move 
upstream  (if  the  surface  wave  moves  downstream).  However,  this  secondary  flow 
is  actually  composed  of  a  superposition  of  exponentially-deoaying  standing 
waves  of  form  ▼Bvo(fi)exp(-0x)sin(Py)exp(-ittt)  where  0  and  »  are  real. 

The  surface  wavinesses  considered  are  either  sinusoidal  waves  of  form 


h(x,t)«h#sin[a(x-ct)]  which  travel  at  speed  c,  or  they  are  galloping  waves  of 
form  htx.tJ^A+Bx+CtlcosIatx-Ujjt)] .  This  wall  shape  induces  a  velocity  of 
form  v(x.t)wvAcosa(x-l^st)  which  propagates  at  the  freestream  speed. 


The  solution  with  wavy  or  galloping  walls  is  superposed  on  the 
corresponding  solutions  for  freestream  disturbances  encountering  a 
semi-inf inite  flat  plate.  The  phase  speed  and  x-wavenumber  of  the  surface 
wave  are  matched  with  the  phase  speed  and  x-wavenumber  of  the  freestream 
disturbance.  Additionally,  the  phase  and  amplitude  of  the  surfaoe  waviness 
are  linked  to  the  phase  and  intensity  of  the  freestream  disturbances.  The 
freestream  disturbance  can  propagate  at  any  speed,  c.  The  result  of 
superimposing  the  wavy-wall  solution  and  the  freeatream  disturbance  solution, 
and  properly  adjusting  phase,  wavenumber,  amplitude,  and  propagation  apeed  of 
the  wavy  wall  is  that  the  family  of  standing  waves  is  eliminated,  along  with 
the  pressure  pulse  at  the  leading  edge.  The  traveling  wave  solution  which 
remains  is  also  modified  from  the  flat  plate  eaae,  and  takes  on  the  form  of 
the  original  freestream  disturbance. 


The  analysis  is  an  example  of  using  traveling  wavy  surfaces  to  modify 
(and  in  this  case  cancel)  standing  waves  and  the  secondary  flow  which  normally 
arise  when  freestream  disturbances  encounter  a  leading  edge. 
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1.  INTRODUCTION 


1.1  Introductory  Cement  g 

Wavy  surfaces  in  engineering  devices  and  nature  do  not  extend  to  infinity 
upstream  and  downstream,  of  course.  They  have  leading  and  trailing  edges, 
they  are  bounded  by  noa-wavy  surfaces,  or  they  connect  with  transition  regions 
which  smoothly  flare  the  wavy  surfaoes  into  rigid,  smooth  surfaces.  A  general 
motivation  for  this  study  is  to  be  able  to  predict  the  effects  of  these  edges, 
boundaries,  or  transition  regions. 

A  second  motivation  is  associated  with  our  awareness  based  on 
boundary-value  problems  by  Tsugd  and  Rogler  (Ref.  1)  that  at  least  five 
families  of  2-D  disturbances  exist  in  a  visoous,  parallel-flow  boundary  layer 
along  a  rigid,  flat  wall.  These  families  include 

(1)  the  discrete  set  of  eigenmodes 

(2)  exponentially-decaying  standing  waves 

(3)  exponentially-growing  standing  waves 

(4)  a  downstream-propagating  continuous  spectrum 

(5)  an  upstream-propagating  continuous  spectrum 

If  the  wall  has  an  active  waw  surface,  where  the  wave  motion  can  be  considered 
to  be  specified  and  the  wavenumber  and  phase  speed  do  not  duplicate  any  of 
those  parameters  of  families  (1)  and  (4),  then  a  wave  (6a)  will  be  found  in 
the  fluid  which  is  associated  with  these  specified  boundary  conditions. 
However,  if  the  wavy  surface  were  compliant  and  two-way  interaction  could 
occur  between  the  materials  on  either  side  of  the  interface,  then  this  forced 
wave  would  become  (6b)  an  eigenwave  that  would  prbpagate  at  nearly  the  speed 
of  the  free  wave  of  the  surface.  These  are  the  Class  B  waves  as  categorized 
by  Benjamin  (Ref.  2).  Also  the  fundamental  eigenwave  from  family  (1)  above 
would  be  modified  by  such  a  compliant  surfaee,  and  this  wave  forms  Class  A 
waves  in  Ref.  2.  If  two  fluids  are  at  different  speeds  and  a  shear  layer 
forms  in  the  interfaoial  region,  then  (8)  the  Kelvin-Helmholtz  instability 
wave  is  a  possible  waveform  and  have  been  called  Class  C  waves.  These 
compliant  or  passive  eases  are  not  further  considered.  The  wall  motions 
considered  in  this  report  are  specified,  although  in  Section  5,  those  motions 
are  linked  to  the  freestream  disturbances. 

In  an  inviscid,  uniform  mean  flow  (U"l) ,  waves  2,3 ,4, and  6  remain  as 
possible  fluctuations  over  an  active  wavy  wall.  The  versions  of  these  waves 
based  on  the  inviscid,  uniform  flow  Rayleigh  equation  can  be  expressed  in 
simple  analytical  forms.  They  have  many  properties  similar  to  those  waves 
found  in  parallel-flow,  viscous  boundary  layers,  although  other  isiportant 
features  of  critical  layers,  sublayers,  and  production  of  unsteady  vorticity 
inside  the  boundary  layer  are  completely  absent  in  the  uniform  mean  flow, 
inviscid  descriptions.  These  simplified  versions  of  the  waves  illustrate  how 
a  wavy  surface  with  a  leading  edge  can  initiate  waves  2  and  6  along  the  plate. 
A  wavy  surface  also  can  modify  the  traveling  wave  solutions  associated  with 
the  freestream  disturbances. 

The  assumptions  of  an  inviscid,  uniform  mean  flow  have  been  introduced  in 
many  analytical  studies  in  the  past.  Irrotational  incompressible  flow  over  a 
small-amplitude  sinusoidal  wall  was  analyzed  by  Kelvin  and  Helmholtz  (Ref.  3). 


This  solution  was  extended  by  Aekeret  (Ref.  4)  to  the  linearized  eases  of 
compressible  subsonic  and  supersonic  flow.  Transonic  flow  over  a  wavy  wall 
was  analyzed  by  Hosokawa  (Ref.  S)  who  included  the  nonlinear  transonic  ten  in 
bis  analysis.  He  assumed  that  tbe  flow  was  isentropic  and  periodic.  Hoore 
and  Gibson  (Ref.  6)  and  Vincenti  (Ref.  7)  analyzed  inviscid  compressible  flow 
over  a  sinusoidal  wall  of  a  gas  in  vibrational  or  chemical  non-equilibrium. 

Horloek  (Ref.  8)  analyzed  flow  in  an  unsteady  wind  tunnel  with 
small-amplitude  sinusoidal  wavy  walls  with  various  phase  relationships  between 
the  two  wavy  walls.  Rogler  (Ref.  9)  analyzed  the  small- amplitude  eases  of 
wall  amplitudes  varying  as  h(x,t)  -  k0ezp[ia(z-ct)]  with  a  complex  and  ctf^, 
and  h(x,t)  ■  (h#'+hxx+hat)exp[ia(x-I^#t)] .  The  latter  wall  ahape  is  a  surface 
wave  which  induces  neutral  disturbances  that  propagate  at  the  freestream 

speed,  c*Uw.  Large-amplitude  traveling  wavy  walls  in  a  channel  were  studied 
in  Ref.  9  also.  The  nonperiodic  flow  induced  by  small-amplitude  traveling 
sinusoidal  walls  A  corner  was  analyzed  by  Rogler  (Ref.  10) . 

The  assumptions  of  an  inviscid,  uniform  mean  flow  have  also  been  usefully 
introduced  in  steady  and  unsteady  linearized  airfoil  theory,  and  the 
propagation  and  dampening  of  sound  in  ducts.  The  uniform  flow  analyses  by 
Benjamin  (Ref.  2)  and  Landahl  (Ref.  11)  were  used  to  study  several  baaic  waves 
in  flows  near  compliant  surfaces. 

A  companion  study  by  Rogler  (Ref.  12)  considers  flow  over  a  semi-infinite 
plate  with  stationary  waviness.  References  associated  with  eight  other 
approaches  which  have  been  used  to  study  flows  over  wavy  walls  are  outlined  in 
that  report. 


l .2  Objectives  aad  pygaaiMti? a  ai  IUa  iaalYtigal/a«gxlffAl  ttafly 
The  objectives  of  this  analytical /numerical  study  are  to 

(a)  Extend  the  previous  analysis  of  flow  over  a  semi-infinite  stationary 

wavy  wall  (Ref.  12)  to  the  case  of  a  semi-infinite  travel ins  wavy 
wall.  This  analysis  also  represents  an  extension  of  the 

Kelvin-Helmholtz  solution  for  flow  over  a  doubly- inf ini te  wavy  wall 
(-••<x<«e)  to  the  semi-infinite  wavy  plate  (0<x<»)  with  a  leading  edge 
at  x«0. 

(b)  Summarize  and  extend  several  classes  of  linearized  boundary 

conditions  of  an  inviscid,  uniform  mean  flow  past  sinusoidal  walls 
with  wall  shapes  of  form  h(x, t)Bhesin[u(x-ct)] ,  and  flow  past 
"galloping"  walls  with  wall  shapes  of  the  form 
h(x,t)»lA+Bx+Ct]cos[u(x-Dc,t)] . 

(c)  Summarize  and  extend  the  solutions  for  freestream  disturbances 
encountering  a  semi-infinite  plate. 


(d)  Combine  the  following  two  solutions  into  a  single  theory 
e  the  solution  for  flow  with  a  semi-infinite  wavy  wall 


e  the  eolation  for  freestream  disturbances  encountering  a 
seai-inf inite  plate 

The  phase  speeds  and  x-wavenmbers  of  the  two  cases  are  the  sane, 
and  the  amplitudes  and  phases  are  adjusted  so  that  the  dipole-like 
pressure  pulse  at  the  leading  edge  is  eliminated. 

The  analysis  serves  as  an  example  of  a  theory  vhich  combines  the  effects 
of  freestream  disturbances  and  surface  vaviness  into  a  single  snalysis.  The 
purpose  of  the  traveling  surfaoe  wave  is  to  modify  or  eliminate  some  of  the 
precursors  of  the  Tollmien-Sehlichting  stability  vaves.  Rather  than  setting 
up  eigenvalue  problems  vhich  consider  the  effects  of  vaviness  on  stability 
vaves*  this  analysis  considers  the  effect  of  vaviness  on  tvo  other  families  of 
solutions:  (1)  the  exponentially-varying  standing  vaves  (represented  here  as 
the  counterpart  described  by  the  Rayleigh  equation  for  a  uniform  mean  flow) 
and  (2)  the  downstream-travel  ing  continuous  spectrum.  No  stability  vaves 
appear  in  this  problem  vith  a  uniform  mean  flow. 

The  olan  of  this  report  is  that  linearized  wall  boundary  conditions  for 
wavy  walls  and  galloping  vails  are  outlined  in  Section  2.  The  solution  for 
flow  past  semi-infinite  plates  vith  wavy  and  galloping  surfaoes  is  presented 
in  Section  3.  The  solution  for  flow  past  a  semi-infinite  £lh£.  plate,  vith 
disturbances  in  the  freestream  (including  arrays  of  vortices  and  irrotational 
fluctuations)  is  summarized  in  Section  4.  The  case  of  an  oscillating  vortex 
sheet  convecting  downstream  past  the  leading  edge  is  also  analyzed  in  Section 
4.  The  combined  case  of  a  wavy  surface  and  freestream  disturbances  is 
analyzed  in  Section  5  by  superimposing  the  solutions  from  Sections  3  and  4. 
The  results  are  discussed  in  Section  6. 


2.  LINEARIZED  WALL  BOUNDARY  CONDITIONS  FOR  FLOW  PAST  WAVY  AND  GALLOPING  WALLS 


In  the  first  part  of  this  section,  s  direct  analysis  will  be  nsed  to 
recower  the  usoal  linearized  boundary  condition  for  an  irrotational  flow  past 
a  small-amplitude  sinnsoidal  wall.  In  the  second  part,  the  case  of 
disturbances  traveling  at  the  freestream  speed,  c~UW(  is  analyzed  by  an 
indirect  method.  The  more  generalized  wall  shapes  which  result  are  called 
"galloping". 


2.1  Sinusoidal  wavy  walls  o£  form  h(z.t)  -  h0ain[a(x-ct)l 

The  geometry  and  notation  for  the  semi-infinite  plate  with  a  wavy  surface 
is  shown  in  Figure  1.  The  displacement  of  the  wall  is  the  traveling  sine  wave 

h(*,t)  -  Wn[a(x-ct]  (2.1) 


where  ha  is  the  amplitude  of  the  surface  waviness,  z  is  the  coordinate  in  the 
streamwise  direction,  a  is  the  x-wavenumber ,  c  is  the  phase  speed,  and  t  is 
time.  The  relation  between  the  surface  displacement  h  and  the  normal  velocity 
v'**  along  the  x— axis  induced  by  that  wall  displacement  and  motion  can  be 
expressed  by  the  linearized  relation 

v^  (x,y*0,t)  ■  dh/dt  +  U^aii/dx  ■  4^^  (2.2a) 

For  the  wall  motion  given  by  eqn.(2.1),  this  velocity  is 

v^(x,0,t)  “  vJ^costaU-ct)]  (for  x>0)  (2.2b) 

where  the  velocity  amplitude  and  wavy  wall  amplitude  are  related  by 


v£^  ■  ah.dJjo-c) 


(2.3a) 


z  «*  x  +  iy 


Figure  1.  Geometry  and  notation  for  flow  past  s  wavy  or  galloping 
semi-infinite  plate. 


terns  of  the  streamf  unction,  this  boundary  condition  is 


^^(x2o,0)  “  h« (I^#-c>sin[a(x-ct)]  +  fnnction(t)  (2.3b) 

is  boundary  condition  along  the  x-axis  downstream  of  the  leading  edge  will 
used  in  Sections  3  and  4  for  sinnsoidal  waves  which  propagate  at  speed  c. 

Although  the  boundary  conditions  (2.3a,b)  are  valid  for  c*!^,  sinusoidal 
rface  waves  in  a  uniform  mean  flow  which  propagate  at  the  freestream  speed 
not  induce  any  flowfield,  y(i)«o.  More  generalized  wall  shapes  which 
duce  disturbances  that  propagate  at  speed  c-U^  are  considered  in  the 
llowing  paragraphs. 


2  Galloping  surfaces  of  fora  h(x,t)«=  [A+Bx+CtJcostaCx-U^t)]  which  induce 
ill propagating  «_£.  speed  c  »  UM 

Some  surface  wave  of  unknown  geometry,  h(x,t),  will  be  assumed  to  induce 
e  velocity  along  the  x-axis 

v(i)  »  v^i^cos[a(x-G00t)]  (2.4) 

sed  on  the  linearized  relation  (2.2a)  between  displacement  and  velocity,  the 
[nation 


dh/dt  +  Uwah/dx  =  v0oos  [a(x-D00t]  (2.5) 

t  a  forced  equation  which  will  admit  both  homogeneous  and  particular 
>lutions.  Only  the  particular  solution  was  considered  in  Section  2.1. 

The  homogeneous  solution  is 

h(H) (x.t)  -  F(x  -  U^t)  (2.6) 

tere  F  is  any  continous  function  with  a  continous  first  derivative  with 
ifficiently  small  amplitude  and  slope  such  that  the  linearized  relation  (2.5) 
t  valid.  This  wall  shape  propagates  to  the  right  at  speed  without 

tanging  shape.  It  does  not  induce  any  disturbances.  This  geometry  can  be 
turier  analyzed  and  represented  as  a  series  or  integral  of  waves  of  forms 

cos[y (x-U^t)]  and  sin[y(x-U#pt)]  (2.7) 

tere  y  is  the  x-wavenumber.  Alternately,  these  are  two  examples  of  the 
mtogeneous  solutions. 

By  the  method  of  undetermined  coefficients,  the  particular  integral  could 
i  assumed  to  be  the  same  form  as  the  forcing  function  and  its  derivative 

Acos[a(x  -  U^t)]  +  Dsin[a(x  -  Uwt)] 

jwever,  these  solution  forms  duplicate  those  of  the  homogeneous  solution 
1.7)  when  Fourier  analyzed  and  y  *  u,  and  must  be  multiplied  by  powers  of  x 
id  t.  The  trial  solution  for  the  particular  integral  is  therefore 


«  [Bx+Ct Jcos [ofx-D^t) ]  +  [Ex+Ft ] s in [a (x-IJ^t ) ] 
By  substituting  eqn.(2.8)  into  eqn.(2.5),  then  the  constants 


(2.8) 


(2.9a) 


and  tbe  two  remaining  constants  are  related  to  the  amplitude  vj1'  by  the 
relation 


n^B  +  C  -  v^i) 


(2.9b) 


in  agreement  with  eqn.(9.11)  of  Ref.  9.  For  a  given  velocity  v£*)f  either 
constant  B  or  C  can  be  specified  and  the  other  constant  found.  B  or  C  can  be 


zero. 


The  complete  solution  is  the  sum  of  the  homogeneous  solution  (2.6)  and 
the  particular  integral  (2. 8, 2. 9a, b) 


h(x,t)  “  Ftx-U^t)  +  [Bx+Ctlcoslafx-Ujjt)] 
The  homogeneous  solution  can  be  rewritten  as 

F(x-D00t)«G(x-D#0t)  +  Acos[a(x-U00t)] 


(2.10) 


(2.11) 


where  A  is  the  appropriate  Fourier  coefficient;  the  function  G(x-Dgot)  is  what 
remains  from  the  homogeneous  solution  after  the  wave  Acosa(x-Uwt)  is 
subtracted  from  F.  In  this  rearranged  form,  the  complete  solution  is 


h(x,t)  *  G(x-U0#t)  +  [A+Bx+Ct]costa(x-l^t)] 


(2.12) 


The  role  of  the  constant  A,  which  is  part  of  the  homogeneous  solution,  is  to 
provide  control  over  the  initial  amplitudes  of  the  surface  waves. 

Neglecting  the  function  G,  several  cases  of  solution  (2.12)  are  possible. 
They  include  the  case  00,  which  represents  waves  decaying  or  growing  in  a 
linear  envelope  in  space,  and  B-0  which  represents  waves  decaying  or  growing 
in  time,  also  in  a  linear  fashion.  The  general  case  has  both  B  and  C  nonzero. 
The  surface  geometries  for  these  galloping  cases  are  plotted  in  Figures  2-4 
from  Ref.  9.  These  figures  help  illustrate  the  role  of  the  constant  A.  Only 
for  some  interval  of  x  (if  00)  or  only  for  some  interval  of  time  (if  B**0) 
will  the  amplitudes  be  sufficiently  small  such  that  the  linearized  boundary 
condition  is  applicable.  The  constant  A  controls  that  interval. 


5.  COMBINED  EFFECT  OF  SURFACE  WAVINESS  AND  FREES TREAM  DISTURBANCES 

ON  A  SEMI- INFINITE  PLATE 


The  linear  solutions  presented  earlier  can  be  superposed.  The  cases  of 
face  waviness  and  the  freestream  disturbances  both  have  the  dimensionless 
s  of  the  normal  velocity 

v(x,0,t)  ■  cos[n(x-t)l 


both  have  the  same  phase.  The  dimensional  x-wavenumber  and  the 
ensional  phase  speed  also  are  required  to  be  the  same  for  the  two  cases, 
ally,  the  characteristic  amplitudes,  v[&)  and  v«^  must  be  the  same  for  the 
cases  so  that  the  characteristic  value  of  the  streamfunctions  are  the 
e.  Dimensionally,  this  requires  that 

-  ohod^-c)  -  v$a)  (5.1) 

n  eqn.(2.3a)  if  the  wall  is  sinusoidal  and  c  and 

v*1*  -  U^B  +  C  -  v<a)  (5.2) 

m  eqn.(2.9b)  if  the  freestream  disturbance  propagates  at  dimensional  speed 
“  U^.  The  constants  B  and  C  are  described  in  Section  2.2;  they  control  the 
elopes  of  the  galloping  wall. 


If  the  above  conditions  are  satisfied,  then  the  streamfunction  for  the 
ibined  effect  of  freestream  disturbances  and  a  wavy  or  galloping  snrface  is 


^(combined)  m  ^(f)  +  p(i) 

m  [jtf(*)  -  ft(i)]  +  fg(i) 

^(combined)  .  ^(a) 


(5.3) 


>  with  this  strategy  of  introducing  a  wavy  wall  or  a  galloping  wall  when 
lestream  disturbances  are  present,  £2  adjustment  is  necessary  for  the 

la.  satisfy  impermeability  At  £he  mt£A£S.  2±  Ol£.  Plate. 

This  result  can  also  be  found  directly  by  consideration  of  the 
indary-value  problem,  without  necessity  to  find  the  solution  (3.3). 
ever,  in  any  practical  device,  the  conditions  on  phase,  x-wavenumber,  phase 
ed,  and  amplitude  cannot  be  exactly  satisfied.  Calculations  with 
matches  are  possible  by  using  the  theories  presented. 

In  summary,  if  the  phase,  x-wavenumber,  phase  speed,  and  amplitude  of  the 
face  waviness  are  adjusted  as  indicated  relative  to  the  freestream 
iturbance,  and  if  the  other  conditions  for  the  theories  for  the  effects  of 
face  waviness  and  freestream  disturbances  are  valid,  then  the  freestream 
iturbance  will  propagate  past  the  plate  as  if  the  plate  were  removed  from 
i  flowfield.  The  iaq>ortant  result  is  that  the  dipole-like  pressure  pulse  at 
i  leading  edge  is  eliminated.  The  theories  provide  the  tools  useful  to 

ermine  the  effects  of  not  exactly  satisfying  the  necessary  conditions. 


igure  8.  Disturbance  streamlines  for  an  oscillating  vortex  sheet  convecting 
downstream  and  alongside  a  semi-infinite  flat  plate.  Three  values  of 
are  used  in  this  figure.  Details  of  the  flow  beneath  the  plate 
and  downstream  of  the  leading  edge  are  plotted  upside  down  in  Fig.  7 
for  t  *  0.5. 


As  noted  in  eqn.(4.2),  the  streamf unction  for  this  vortex  sheet 
teracting  with  s  semi- inf inite  piste  is  given  by  -  ^(s)  -  ^(i),  where 

is  given  by  eqns. (4.4b #4 .5b) .  The  resultant  streamline  pattern  is 
otted  in  Figure  8. 

In  the  next  section,  the  cases  of  the  flow  induced  by  surface  waviness 
id  the  flow  associated  with  freestream  disturbances  interacting  with  a 
imi-inf inite  plate  are  superposed. 
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•  Stationary  surface  waviness  (Kef.  12) 

e  Traveling  surface  vaviness  (Sections  2 >3  of  this  report) 
e  Combination  of  traveling  surface  vaviness  and  freestream  disturbances 
(Section  5  of  this  report) 

e  Quiescent  freestreaa  with  a  flat  plate  (v^wv^^-O) . 


A 


4.2  Interaction  si  1A  oaolllatins  vortex  thC.C.t  llti  A 


aJLaAS. 


As  an  exaaple  of  the  procedure  suaaarised  above,  the  problem  of  an 
oscillating  vortex  sheet  lying  along  the  line  y  ■  y,  will  be  analyzed.  This 
sheet  oonveots  downstream  at  speed  UH  past  the  leading  edge  of  the 
semi-infinite  plate. 

In  order  to  define  the  characteristic  velocity  induced  by  the  vortex 
sheet,  the  equations  of  this  section  are  in  dimensional  form.  If  a  vortex 
sheet  of  strength  r  per  unit  length 


y(x,t)  -  r«sin[a(x-t)l 


(4.3) 


eonveets  along  the  line  y  *  y, ,  then  the  streamfunction  induced  by  that  vortex 
sheet  is 


j^(a)  m  C1e“oysin[o(x-t)]  (for  y  >,  y,) 
/<•>  -  C,e+«7sin[a(x-t)]  (for  y  <  yt) 


(4.4a) 


(4.4b) 


where  C,  and  C,  are  constants  which  will  be  linked  to  the  amplitude  of  the 
vortex  sheet,  ra. 

Since  must  be  continuous  across  the  sheet,  then 


C,exp (-ay, ) sin [a (x-t ) ] *C,exp (ay* ) sin [a (x-t ) ] 
or  Cx  ■  C,exp(2ay,) 

Across  the  vortex  sheet,  the  jump  condition  must  be  satisfied 
[_^(a)]  _  r_^a)i  _ 

y  Jr,+  1  *  yJy,-  T 

or  C,  m  r, exp (-ay,) 


(4.5a) 


(4.5b) 


If  the  vortex  sheet  lies  above  the  plate,  then  y,  >  0,  and  the  normal  velocity 
along  the  x-axis  is 


v(*)  ■  C,ocos [a(x-t) ] 


or  the  characteristic  velooity  is 


▼e*)  "  <*T#»*P<-«ys) 


(4.6) 


This  ampl itude  is  required  in  the  nondimens i one 1 ization,  as  noted  in  the 
sentence  following  eqn. (4.2). 
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4.  INTERACTION  OF  FREESTXEAM  DISTURBANCES  WITH  A  SEW- INFINITE  FLAT  PLATE 


r 

i 


4.1  Introduction 

This  section  summarises  the  analyses  of  the  inviscid  interaction  of  2-D 
freestream  disturbances  with  a  son i- inf inite  flat  plate.  Hie  disturbances  can 
have  rotational  and/or  irrotational  regions.  Different  disturbances  propagate 
downstream,  upstream,  or  are  stationary  in  tine.  A  number  of  studies  have 
been  carried  out  wbieb  are  listed  below: 

e  Array  of  square  vortices  oonvecting  downstream  (Ref.  14) 
e  Array  of  rectangular  vortices  couvecting  downstream  (Ref.  IS) 
e  Oblique  and  longitudinal  plane  waves  of  vortieity  oonvecting  downstream 
(Ref.  IS) 

e  Oscillating  vortex  sheet  oonvecting  downstresm  (Section  4.2  of 
this  report) 

e  One-sided  irrotational  wave  propagating  at  speed  c  (Ref.  IS) 
e  Double-sided  irrotational  waves  with  oblique  or  normal  velocities 
along  the  x-axis  and  propagating  at  speed  c  (Ref.  IS) 
e  Double-sided  irrotational  wave  whioh  produoes  only  longitudinal 
velocities  along  the  x-axis  (Ref.  IS,  Note  that  v(») (x,0,t)-0) 
e  Row  of  potential  vortices  (Ref.  IS) 

e  Karman  vortex  street  propagating  at  speed  c  downstream  (Unpublished) 
e  Half-plane  of  rectangular  vortices  couvecting  downstresm  (Unpublished) 

With  the  exception  of  the  row  of  potential  vortices  and  the  Karman 
street,  the  above  disturbances  have  the  velocity 

v**)(x.O,t)  ■  vJ4^cos[a(x-ct)]  (along  the  x-axis)  (4.1) 

The  superscript  (a)  emphasizes  that  this  is  the  velooity  of  the  freestream 
disturbance  in  the  absence  of  the  plate,  i.e.  if  the  plate  were  removed  from 
the  flowfield.  What  is  the  flowfield  with  the  plate  present? 

Under  conditions  discussed  in  Refs.  14,15,  if  {A»)  is  the  streamf unction 
with  the  plate  absent,  then  the  streamf unction  with  the  plate  present  is 

j4(f)  .  $*(a)  _  jt(i)  (4.2) 

where  (£(*)  it  the  streamf unction  given  in  eqn.  (3.3).  The 

nondimensional izations  are  the  same  as  before,  except  that  the  characteristic 
velooity  is  the  amplitude  of  the  normal  velocity  along  the  x-axis  of  the 
disturbance  in  the  absence  of  the  plate,  v£b)(  lt  note(j  in  8qn.(4.1). 

The  Karman  street  and  the  row  of  potential  vortices  traveling  downstream 
have  velocities  whioh,  in  the  absence  of  the  plate,  are  periodio  in  x,  but 

they  are  not  sinusoidal.  The  interactions  of  these  disturbances  with  s 

semi-infinite  plate  have  been  expressed  as  Fourier  series. 

Hence,  the  streamfunotion  {AD  plays  a  role  in  the  interaction  of  a 
variety  of  freestream  disturbances  with  a  plate,  as  well  as  the  flowfield 
induced  by  wavy  and  galloping  surfaces.  To  the  above  list,  the  following 

cases  oan  be  added 


v(,)(x>0,y,t)  -  ▼$*)<p)«ia0y«xp(-px-it»t) 


(3.8) 


i-y.i 


where  p  is  «  real  mtnnbii.  Each  of  these  waves  satisfies  iaperaeability, 
each  decays  exponentially  in  the  streaawise  direction,  and  each  oscillates 
sinusoidally  in  tiae.  A  difference  between  the  the  waveform  (3.8)  above  and 
the  results  of  Kef.  12  is  that  the  above  flow  is  unsteady,  while  Ref.  12  is 
steady. 
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An  analogous  result  for  freestreaa  disturbances  encountering  a  plate  was 
presented  in  Ref.  Id.  An  effect  of  the  leading  edge  was  to  excite  a  speotrna 
of  standing  waves,  and  also  aodify  the  trsveling  wave  solution  for  the 
freestreaa  disturbance. 


In  conclusion  of  this  section,  the  stresafunetion  and  velocities  of  the 
flow  induced  by  a  seal- inf inite  wavy  or  galloping  wall  have  been  found  which 
illustrate  the  flow  around  the  leading  edge  and  other  upstream  influences. 
Downstreaa  of  the  leading  edge,  the  flow  is  represented  in  terns  of  the 
classical  Eelvin-Helaholtz  solution  and  a  secondary  flow.  This  unsteady, 
irrotational,  secondary  flow  satisfies  the  taperne ability  condition  on  a  flat 
plate,  and  can  be  represented  as  a  superposition  of  exponentially-decaying 
standing  waves. 
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f(x)  m  -  S(z)]oo«(jtz  /2)  -  [-J —  C(z)]sin(nz*/2) 

g(z)  ■  C(z)]cos(flz*/2)  +  t-J--  S(z)]iin(nz*/2) 

S(z)  -  Sa(nz*/2) ;  C(z)  - 

The  asymptotic  series  for  f  and  g  (Kef.  12,  page  302)  are 


itzf(z)~  1  +  £(-l)»  [1.3...(4sf-l)/(jtz*)*a] 

m 

nzg(z)~  £(-l)“  I1-3...<4b+1)/(rz*)4“+1] 

These  series  converge  rapidly.  The  solution  form  (3.5d)  avoids  the  problems 
of  finding  the  small  differences  between  large  numbers.  While  the  direct 
evaluation  of  solution  (3.5b)  was  limited  to  about  |z|*1.8  when  32  bit 
arithmetic  was  used,  and  was  limited  to  about  |z|-5.2  when  64  bit  arithmetic 
was  used,  solution  (3.5d)  has  been  evaluated  at  |z|*400,  which  is  200 
wsvelengths  from  the  leading  edge. 

A  singularity  occurs  at  the  leading  edge  for  the  velocities,  with  the 
validity  of  the  theory  consistent  with  linearized  airfoil  theory.  The  region 
where  the  velocities  are  large  oan  be  made  as  small  as  desired  by  making  the 
amplitude  of  the  wall  waviness  smaller.  Henee.  this  is  the  proper  first-order 
solution  as  long  as  the  flow  is  unseparated  at  the  leading  edge.  Sinoe  the 
practical  case  of  flows  past  bodies  with  streamlined  leading  edges  is  of 
interest,  the  assumption  of  unseparated  flow  isqtlieitly  contained  in  Laplace's 
equation  appropriately  models  that  feature  without  unduly  complicating  the 
analysis  with  details  of  the  nose  geometry  and  plate  thickness. 

Solution  (3.3)  with  streamlines  plotted  in  Figs. 6a. b  is  valid  both 
upstream  and  downstream  of  the  leading  edge.  As  noted  in  Kef.  12.  the  leading 
edge  modifies  the  Kelvin— Helmholtz  solution,  denoted  by  ff(kh) .  Hence,  the 
Kelvin-Helsdioltz  solution  is  subtracted  from  (6^ .  and  the  result  is  a 
streamf unction  for  a  secondary  flow 

^(s)  „^(i)  _  j*(kh)  (for  x>o)  (3.6) 

The  streamlines  of  this  secondary  flow  are  shown  in  Fig.  7.  Since  both 
and  zre  irrotational ,  then  is  also  irrotational.  Sinoe  both 

and  have  the  same  wall  boundary  conditions,  then 

(x>0,y-0,t)  “0  (3.7) 

This  result  denotes  that  this  secondary  flow  satisfies  the  ia^ermeabil ity 
condition  along  a  flat  plate.  However,  it  has  indirectly  been  excited  by  the 
combined  effects  of  the  wavy  wall  and  the  leading  edge. 

Based  on  the  mathematical  developments  of  Ref.  12,  the  secondary  flow  can 
be  represented  as  a  superposition  of  exponentially-decaying  standing  waves  of 
the  form 
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Figure  6 


.  Disturbance  streamlines  for  the  flow  induced  by  a  wavy  or  galloping 
surface  of  a  semi-infinite  plate  for  two  different  times  (From 
Ref.  12  except  that  the  phase  angle  of  the  stationary  wariness  is 
replaced  by  irt) . 


13 


"  .  *  4  •  >•»  •  •  »  *  -  *  «  •  •  *  .  *  .  •  . 


>  1  % 


problem  arose  associated  with  the  interaction  of  freeatreaa  disturbances  with 
a  semi- inf inite  flat  plate. 

The  solution  for  the  streamf unction  of  the  flow  induced  by  the 
semi-infinite  plate  with  a  sinusoidal  wall  is 

^(i)  «  {cosnt  Real [-cosnzSa (-nz) -sinnzCi (-nz) 

+cosnzCa  (— jtz)  ‘"sinirzSj  (-nz)  +sinnz] 

-slant  Real [-eosnzSa (-nz) -sinnzC, (-nz) 

-coanzC,  (-nz)  +sinjtzS,  (-nz)  +cosnz] }  /n  (3 .3 ) 

where  Sa  and  Ca  are  the  Fresnel  integrals 


Sa  -  (2n)"l/aJt"l/*«in(t)  dt 


-a  ■  (2n)-1^*Ji 


-s/sf -l/S 


cos(t)  dt 


(3.4a) 


(3.4b) 


The  Fresnel  integrals  are  evaluated  numerically*  either  as  asymptotic  series 
valid  for  large  argument  or  series  expansions  valid  for  small  argument 
(Ref.  13).  The  details  of  the  evaluations  are  given  in  Ref.  12. 

Figures  6a, b  show  the  disturbance  streamline  patterns  for  the  flow 
induced  by  the  semi-infinite  wavy  wail  at  two  times  a  quarter  period  apart. 
These  figures,  and  the  others,  do  not  include  the  uniform  mean  flow.  The 
induced  flow  depends  on  the  time.  In  all  cases,  the  wavy  surfaces  on  opposite 
sides  of  the  plate  are  in  phase.  The  upstream  influence  of  the  waviness  can 
be  clearly  seen,  as  expected  since  Laplaoe's  equation  is  elliptic. 


The  velocities  corresponding  to  solution  (3.3)  are 
n(i)  .  «  -iF(z,t) 

v(i)  .  0(i)  .  F(z,t) 


(3.5a) 


(3.5b) 


where  F  is  the  complex  function 

F  «  sin(nz-nt) [Sa(-nz)-Ca(-nz)l 

-cos  (nz-nt )  [Sa  (-«z)  +C8  (— nz )  — 1]  -( sinnt+cosirt )  / [n (-2z) 1 )  (3.5c) 

For  large  |z|,  computational  experience  has  shown  that  it  is  better  to  recast 
the  solution  as 

F  ■  sinnttf  -  g  -  l/[r(-2z)1^*]}  +  cosntCf  +  g  -  1/ [n(-2z) ) }  (3.5d) 
where  f(z)  and  g(z)  are  related  to  the  Fresnel  integrals  by 
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Figure  5 


.  Conformal  mapping  of  the  upper  half-plane  onto  a  quarter— pi 
the  resnltant  boundary— value  problem. 


3 .  FLOW  OVER  A  SEMI- INFINITE  PLATE  WITH  TRAVELING  WAVES 


The  analysis  of  this  section  is  closely  related  to  the  analysis  by  Rogler 
(Ref.  12)  for  flow  over  a  stationary  wavy  wall,  except  in  the  present  ease, 
the  surface  wave  travels  (downstreaa  or  upstream)  at  speed  c.  The  case  of 
galloping  walls  is  also  included. 

By  conformal  mapping  and  integral  transforms,  the  flowfield  with 
traveling  surface  waviness  or  galloping  on  a  semi-infinite  plate  will  be 
analyzed.  A  solution  is  sought  which  satisfies  Laplace's  equation,  V* 

■  0,  where  is  the  streamf unction  for  the  flow  induced  by  the  surface 
waviness.  For  sign  convention,  the  streamfunction  is  related  to  the 
longitudinal  and  normal  velocities  by  and  respectively. 

The  streamfunction  along  the  z-axis  downstream  of  leading  edge  is 

$£^(x,0,t)  ■  (v£D/a)cos[a(x-ct)l  (for  x>0)  (3.1) 

When  written  in  this  form,  the  theory  is  applicable  for  the  two  cases 
described  in  the  previous  section:  (1)  v£i'  «  ah0UW~c)  for  wavy  walls  of 
sinusoidal  form,  and  (2)  v*1*  .  +  C  for  galloping  walls  which  induce 
disturbances  traveling  at  the  freestream  speed. 

Hereafter,  the  streamfunction  is  nondimensionalized  against  v^n/a  and 
the  velocities  are  nondimensionalized  against  v£i).  The  x  and  y  coordinates 
are  nondimensionalized  against  the  half-wavelength  A/2*n/a  of  the  sinusoidal 
disturbance.  For  historical  reasons,  the  time  is  nondimensionalized  against 
the  characteristic  time  n/ac,  which  is  the  time  for  a  fluid  particle  traveling 
at  speed  c  to  travel  a  half-wavelength,  Xj 2«n/a.  Under  these 
nondimensionalizations,  the  amplitude  of  the  streamfunction  is  1/n,  and  the 
amplitude  of  the  normal  velocity  along  the  plate  is  unity. 

For  the  case  of  the  waviness  being  in-phase  on  the  opposite  sides  of  the 
plate,  the  streamfunction  will  be  symmetric  about  the  x-axis  and  the 
homogeneous  Neumann  boundary  condition  is  applicable  upstream  of  the  plate 

-  0  (for  y-0,  x<0)  (3.2) 

The  solution  is  asstmiod  to  vanish  far-away  from  the  plate,  both  far-upstream 
and  at  large  distances  laterally  from  the  plate.  Far-downstream  of  the 
leading  edge,  the  solution  is  bounded.  It  will  be  shown  that  the  classical 
Kelvin-Helaiholtz  solution  is  recovered.  These  exterior  boundary  conditions, 
along  with  the  Dirichlet  and  Neumann  conditions  along  the  x-axis  and  Laplace's 
equation  yield  a  well-posed  elliptic  mathematical  system. 

The  half-plane  above  the  x-axis  is  conformally  mapped  onto  the 
quarter-plane  as  shown  in  Figure  5  by  the  mapping,  p  “  z1/*,  where  z  *  x+iy 
and  p  “  t+it).  The  conformal  mapping  alleviates  the  problem  of  split  boundary 
conditions  along  the  x-axis,  with  a  Neumann  boundary  condition  for  x<0 
(eqn.  3.2)  and  a  Dirichlet  boundary  condition  for  x>0  (eqn.  3.1).  Under  this 
mapping,  the  resultant  mathematical  system  is  defined  in  Figure  2.  The 
solution  is  obtained  by  a  cosine  integral  transform  in  the  (  direction.  The 
details  are  presented  in  Ref.  14,  where  a  closely  related  boundary-value 


6.  SUMMARY,  DISCUSSION,  AND  CONCLUSIONS 


Analyses  of  linearized  perturbations  from  a  uniform  mean  flow  have  been 
carried  out  for  three  oases: 

(1)  Flow  past  sinusoidal,  traveling  surface  waves  on  a  semi-infinite 
plate,  and  flow  past  more  generalized  "galloping"  wall  shapes  which 
induce  a  neutral  normal  velocity  along  the  plate  which  propagates  at 
the  freestream  speed,  c“Ufl0. 

(2)  Flow  with  freestream  disturbances  encountering  a  semi-infinite 

plate.  An  example  analysis  of  an  oscillating  vortex  sheet 

convecting  downstream  is  given,  and  references  for  many  other 
vortical,  irrotational,  and  mixed  rotational/ irrotational 
disturbances  are  cited. 

(3)  Flow  with  surface  waviness  and  with  freestream  disturbances  is 
analyzed  by  adjusting  the  phase,  x-wavenumber,  phase  speed,  and 
amplitude  (?)  of  the  surface  waviness  or  galloping. 

The  results  for  the  first  two  cases  clearly  illustrate  the  upstream 
influence  of  the  semi-infinite  plate  in  this  elliptic  problem.  The  classical 
irrotational  Kelvin-Helmholtz  solution  is  recovered  downstream  of  the  leading 
edge,  but  an  effect  of  the  leading  edge  is  to  initiate  a  secondary  flow.  This 
secondary  flow  satisfies  the  impermeability  condition  on  a  flat  plate.  This 
secondary  flow  is  also  irrotational. 

This  secondary  flow  appears  in  the  first  two  cases.  However,  for  the 
superposed  cases  of  freestream  disturbances  and  wavy/galloping  surfaces,  the 
secondary  flow  is  eliminated  if  the  phase,  x-wavenumber,  phase  speed,  and 
amplitude (s)  of  the  wavy  or  galloping  surfaee  are  adjusted  to  those  properties 
of  the  freestream  disturbances. 

As  shown  in  Refs.  12,1 6,  this  secondary  flow  is  composed  of  a 
superposition  of  exponentially-decaying  standing  waves  downstream  of  the 
leading  edge.  Hence,  the  analyses  of  semi-infinite  wavy  walls  show  that  the 
motion  of  the  wall  can  excite  this  family  of  waves  if  a  leading  edge  is 
present.  Upstream  of  the  leading  edge,  the  flowfield  around  the  leading  edge 
induced  by  wall  motion  can  be  represented  as  a  superposition  of  growing 
standing  waves. 

This  analysis  is  based  on  linearized  versions  of  the  impermeability 
boundary  condition  at  the  wavy  surfaoe,  and  the  mean  flow  is  assumed  to  be 
uniform.  Just  as  the  Kelvin-Helmholtz  solution,  linearized  steady  and 
unsteady  airfoil  theories,  and  theories  of  uniform  flows  past  (doubly 
infinite)  compliant  surfaces  have  had  a  useful  role  in  the  subsequent 
development  of  those  subjects,  the  introduction  of  these  assumptions  in  a 
first  analysis  is  justified.  Of  course,  such  theories  cannot  account  for 
effects  in  the  boundary  layer,  but  they  may  provide  useful  edge  boundary 
conditions  or  solution  forms  for  the  boundary  layer  analyses.  An  example  of 
this  information  is  the  recognition  that  a  secondary  flow  arises  as  a  result 
of  the  leading  edge. 


Vail  motions  that  induce  disturbance  a  which  propagate  at  the  freestream 
speed  were  found  analytically  in  Section  2,  as  essentially  unchanged  from 
Ref.  9.  These  analyses  were  carried  out  under  the  context  of  a  uniform  mean 
flow,  where  the  usual  sinusoidal  wall  motion  traveling  at  speed  Uw  voui£  aot 
induce  any  disturbances  since  -  0  -  0  in  eqn.  (2.3a).  The  reason  that  the 
case  of  c  a  Uw  is  a  useful  case  is  that  low-intensity  vortical  fluctuations  in 
a  viscous  flow  propagate  at  nearly  that  speed. 

However,  if  a  parallel-flow  boundary  layer  were  present  and  viscous 
effects  were  incorporated,  then  the  linearised  wall  boundary  conditions  along 
the  x-axis  would  be  (e.g.  Ref.  18) 

u(x,0,t)  -  -Oy(0)h# 

v(x,0,t)  -  ch# 

Hence,  the  effect  of  the  mean  boundary  layer  is  to  provide  relative  motion 
between  the  mean  flow  and  the  surface  wave  traveling  at  Uw.  these  boundary 
conditions  indicate  that  disturbances  are  induced  by  a  sinusoidal  wave 
propagating  at  speed  Uw  if  a  mean  boundary  layer  is  present.  Further 
examination  of  the  boundary  conditions  fox  c  near  necessary  to 
understand  the  relationship  between  the  uniform  flow,  inviscid  boundary 
condition  and  the  viscous  boundary  conditions  with  a  mean  boundary  layer  as 
the  Reynolds  number  becomes  infinite  and  all  visoons  layers  have  a  vanishing 
thickness. 

In  conclusion,  the  analyses  for  semi-infinite  wavy  walls  show  that  the 
motion  of  the  surface  wave  induces  a  secondary  flow  near  the  leading  edge 
which  is  superposed  on  the  classical  Kelvin-Helmholtz  solution  for  flow  over  a 
traveling  wavy  surface.  This  secondary  flow  is  composed  of  a  superposition  of 
exponentially-decaying  standing  waves.  When  freestream  disturbances  are 
present,  the  corresponding  secondary  flow  associated  with  the  disturbances 
interacting  with  the  leading  edge  can  be  eliminated  if  the  phase, 
x-wavenumber,  phase  speed,  and  amplitude (s)  of  the  wavy  or  galloping  surface 
are  properly  adjusted.  Further  studies  of  these  effects  on  boundary  layers 
are  warranted,  and  additional  studies  of  traveling  or  galloping  surface  waves 
are  needed  when  the  phase  speed  is  exactly  or  nearly  the  freestream  speed. 
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